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Symbols (unless otherwise noted in report) 


a 


acceleration of a point 


a(t) 



A, B, C, D 
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Co> Cl, C2, C3 


a function of time which increases from 0 ton / 2 gradually 

the coordinates of the upper bearing block of the j-th 
actuator in Frs components 
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NASA adaptive algorithm washout parameters 
mathematical mean of a statistical variable 

actuator extensions 

pilot’s sensation error 

a matrix which relates the states of the optimal control 
system with the input to the controlled plant 

reference frame 


f 

f 


Go 


Gc 


g 

H 


specific force = a - g 
sensed specific force 

a constant in the UTIAS otolith sensation model 
UTIAS adaptive algorithm steepest descent parameters 
acceleration due to gravity 

a matrix related to the white noise input in the state-space 
model of a control system 
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Jsub 
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K$ub 


l^sub» Psub 


Lsi 

P.(t) 

P 
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Ri» Rz* Ra* R 
P» Pi» Q> Rd 


E 

S 

s 

Sq s, 
t 

Ts 


To,Ti,T2,T3 

TF 

u,v,w 


identity matrix 
system cost function 

a constant in the UTIAS otolith sensation model 

NASA adaptive algorithm steepest descent parameters 

UTIAS adaptive algorithm washout parameters 

displacement of the i-th hydraulic actuator 

transformation matrix fix>m fiame £,s into frame Eri 

solution of the algebraic Riccati equation 

roll rate component of ® 

pitch rate component of © 

yaw rate component of ® 

weighing matrices in a cost function 

weighing matrices in a cost function 
radius vector 

simulator centroid displacement 
Laplace variable 

slopes of the nonlinear gain polynomial 
time 

transformation matrix from angular velocity to Euler angle rates 
coefficients in the semicircular canal sensation model 
transfer function 

X-, y-, and z-components of velocity Y 
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m: 

X 

XL 

Y 

W(s) 

x,X 


X 

I 

a 

P 

e 


p 

a 




COb 


On 


O 


input to a control system 

input to a standard form optimal control system 

filtered white noise 

white noise 

velocity Y = [u v w]^ 

optimal control open-loop transfer function matrix 

system state vector 

output of a control system 

output of the pilot’s vestibular system 

rotation angle of the simulator platform 
Euler angles 0 = [ <|) 0 <p 
Euler pitch angle component of 0 
Euler roll angle component of 
Euler yaw angle component of 
scalar weight in optimal cost criterion 
cost criterion for optimal algorithm 

time constants in the semicircular and otolith sensation models 

T 

angular velocity about the body frame fii = [ p q r ] 
first-order high-pass filter break frequency 
second-order system imdamped natural frequency 
sensed angular velocity 
second-order system damping ratio 
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Subscripts (main symbol)subscript 

In most cases, subscripts indicate to what the main symbol is related. 
( )a relates to aircraft 

( )ar relates to aricraft rotation 

( )cA relates to aircraft point corresponding to simulator centroid 

( )cG relates to center of gravity of aircraft 

( )d relates to simulator states included in the cost function 

( )dis relates to filtered white noise disturbance model 

( )des relates to desired output 

( )i relates to inertial reference frame 

( )l relates to low pass filtering 

( )n relates to white noise input states 

( )oio relates to otolith model 

( )pA relates to pilot in aircraft 

( )ps relates to pilot in the simulator 

( )r relates to reference model 

( )s relates to simulator 

( )sR relates to simulator rotation 

( )ssc relates to semicircular canals sensation model 

( )sT relates to simulator tilt coordination 

( )v relates to pilot’s vestibular model 

( x,y, or z components 
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Superscripts (main symbol)*”'^'^"'’* 

In most cases, superscripts indicate which frame the main symbol is in 
( )'^ in^rcraft frame FrA 

( )' in inertial frame Fri 

( )^ in simulator frame Fr$ 
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Abstract 


This project was conducted in two phases. In Phase I, four algorithms including 
the classical algorithm, the NASA adaptive algorithm, the UTIAS adaptive algorithm, 
and the optimal algorithm were investigated. The classical algorithm generated results 
with more distortion, more delay and lower magnitude than the results generated by the 
other three algorithms. The classical algorithm is the fastest one. This is of little 
importance since today’s computers are fast and none of the four algorithms will run 
beyond the required time. Therefore, the classical algorithm has no advantage and was 
not considered in Phase 2. 

The two adaptive algorithms are basically similar. The NASA algorithm is well 
tuned with satisfactory performance. The UTIAS adaptive algorithm strives for more 
flexibility, but results show that it does not behave better than the NASA algorithm while 
having more undesirable properties. Some changes were made to the adaptive algorithms 
such as reducing the magnitude of undesirable spikes. 

The optimal algorithm was found to have the potential to behave much better than 
it did in Phase 1. In Phase 2, the optimal algorithm was redesigned. The center of 
simulator rotation was redefined. More terms were involved in the optimal algorithm 
cost function to yield more flexibility in tuning the algorithm. A new design approach 
featuring a Fortran/Matlab/Simulink interactive design was used. Each set of selected 
parameters could be tested in only 30 seconds while the old design approach could 
require as much as 15 minutes. This makes it possible to try hundreds of sets of 
parameters. As a result, the optimal algorithm was well tuned in Phase 2 that also 
incorporated a revised vestibular sensation model. 
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The effect of motion was also discussed. The topic covers what type of motion is 
desirable and what type is undesirable. A new semicircular canal sensation model was 
constructed and justified. The discussion and the new sensation model helped to develop 
the optimal algorithm and to evaluate the motion-base driving algorithms. 

In Phase 2 comparisons were made between the NASA adaptive algorithm and 
the redesigned optimal algorithm. Results showed that the optimal algorithm has some 
advantages over the NASA adaptive algorithm and might be the best among the four 
algorithms involved in this study. At the same time the NASA adaptive algorithm was 
observed to be a very well developed algorithm. 

There were some general problems left unresolved in Phase 1 that required 
solutions. The first problem was that when inputs had large magnitudes, all the 
algorithms tended to drive the simulator beyond its motion limit. A new nonlinear gain 
algorithm was designed in Phase 2, the effect of which was that the simulator would not 
reach its motion limit in most input cases, while the gain in low magnitude input cases 
remained high. The second problem was that when the simulator reached its limit, an 
algorithm was needed to brake the simulator to a full stop and then release the brake at 
some proper time to allow the simulator to follow the output of the washout algorithm 
again. This braking algorithm was developed in Phase 2. 
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1. Introduction 


With continuing improvement in hardware and software, flight simulation plays 
an expanding role in the training of aircraft crews, design of new aircraft, and 
entertainment. This study evaluates a motion-base driving algorithm for a modem six- 
post synergistic aircraft simulator. Three types of algorithms, classical, adaptive, and 
optimal, are evaluated within the scope of this study. Some implemented cueing 
algorithms were first investigated, with some efforts then spent to improve or even re- 
design them. 

The purpose of a motion simulation is to provide task-critical motion and force 
inf ormation (i.e., cues) and any required components of the stressor-induced workload 
increment that would be present in flight. Since a ground-based flight simulator system 
cannot duplicate the motions of an actual aircraft, it becomes necessary to determine the 
best way to utilize its limited capabilities to provide the most necessary and beneficial 
motion cues. It is also critical for the cueing algorithm to avoid any improper motion 
cues since it is commonly known that improper motion cues in some flight conditions 
have great negative effects on the simulation. It is reported that some motion systems 
experienced being turned off to avoid improper motion cues. A principle component of a 
motion simulator design is the determination of the motion information that is relevant to 
the task, can impact human performance, and can be provided within technical and 
economic constraints. This requires some knowledge regarding the human’s motion 
perception that is therefore also an important portion of this study. 

The motion system usually works in conjimction with a visual system to 
accomplish effective simulation. Although some of the aspects of the visual system are 


11 



taken into consideration in this study, a detailed study of the visual system or the co- 
operation of the two systems is beyond the scope of this research. 

This project was conducted in two phases. In Phase 1, four motion cueing 
algorithms, the classical algorithm, the NASA adaptive algorithm, the UTIAS adaptive 
algorithm, and the optimal algorithm were investigated. The results of these algorithms 
were compared with modifications made to correct problems observed with both the 
NASA adaptive algorithm and the UTIAS adaptive algorithm. In Phase 2, the NASA 
adaptive algorithm and a redesigned optimal algorithm were further investigated, with the 
optimal algorithm incorporating a revised model of the human vestibular system. 
General problems common to all algorithms not resolved in Phase 1 were also addressed 
in Phase 2. 

The motion cueing algorithms are intended to drive the MCFADDEN 676B-B046 
simulator at the NASA Langley Research Center. The simulator is a hydraulic six-post 
synergistic full motion simulator. 
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2 . Background Information 

2.1 . Effect of Simulator Motion and Human Motion Perception 

The goal of the motion system, along with the visual system, is to provide a 
vir tual environment to the pilot so that in a simulator the pilot can perform the controls 
and maneuvers which are consistent with how they are to be performed in a real aircraft. 
The simulation can be employed to help in the design of new aircraft, the training of 
pilots, and research. What is critical to understand is whether motion helps the 
simulation and what kind of motion is really desirable. Although early studies showed 
some controversy on the effect of motion on aircraft simulation, results of further studies 
converged to some consistent conclusions. 

2.1.1. Effect of Motion Cue 

Gundry [1] reported that Douvilier, et al., Matheny, et al.. Perry and Naish, and 
Tremblay, et al. foimd that although motion did not always help to reduce pilot 
performance error in some simulation tests, there were always differences in the pilot 
control activity power spectra. When motion was provided, there was an increase in the 
occurrence of high-frequency/low amplitude control movements. These changes served 
to make the control activity in the moving simulator appear more like that observed 
during flight than that recorded in a fixed base simulator. The presence of motion was 
found to reduce phase lag, increase the mid-frequency gain and crossover frequency, and 
reduce the size of the remnant. 

Gundry [1] also reported that Sadoff, et al., and Meiry found that when the 
simulator dynamics were stable, the presence of maneuver motion did little to improve 
control. But as the vehicle became unstable, maneuver motion became more important; 
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its presence allowed the operator to exercise control even in regions where control by 
visual cues alone would be impossible. 

The visual system alone could provide motion illusions in many simulation tests. 
It was foimd that as long as there was some proper simulator motion at the beginning of 
each maneuver, the motion illusions introduced by the visual system could be established 
more easily and faster. In some simulation cases the motion illusions could exist all the 
time after its establishment without continuing involvement of simulation motion. This 
implies that the simulation motion onset requires high attention. 

Improper simulator motion could be very harmful for it might conflict with the 
visual cue and then break the motion illusions introduced by the visual system. The 
motion cue that conflicts with the visual cue is called a negative motion cue in the 
literature. Negative motion cues should always be avoided whenever possible. 

Clark, et al. [2] found that a pilot’s vestibular system could process low level, 
constant acceleration in the presence of vibratory acceleration as efficiently as it could 
without the vibratory noise. In other words, vibratory acceleration had little or no 
masking effect on the detection of constant acceleration over a wide range of intensity 
levels of constant acceleration. This implies that it is not practical to employ vibrations 
to mask some motion sensation such as the rotational sensation when the simulator tilts to 
simulate sustained linear acceleration. 

2.1.2. Human Motion Perception 

As pointed out in the UTIAS report [3], [4], deriving the human’s motion 
sensation models is important for the evaluation of motion-base drive algorithms and the 
formulation of the optimal control algorithm. The human vestibular system located in the 
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head is found to be dominant in human motion sensation. The vestibular system consists 
of two important parts. One part is the semicircular canals that sense rotational motion 
and the other are the otoliths that sense linear motion. The UTIAS report [3], [4] 
provided models for both the semicircular canals and the otoliths. 

The semicircular canal sensation model is given as: 

2 

03 _ 

ro ■ (•C,S+l)(TsS+l)(T.S+l) 


where is the angular velocity input and ra is the sensed angular velocity. Tl, Ts, and Xa 
are time constants. The model represents a second order torsion pendulum with an 
adaptation term. Tl has a unique value for each rotational degree of freedom: 6.1 sec (roll 
input p), 5.3 sec (pitch input q), and 10.2 sec (yaw input r). t$ = 0.1 sec and Xa = 30 sec 
for all three degrees of freedom. 

The otolith sensation model is given as: 


f_ K(x,s + 1) 

f (XlS + 1)(XsS + 1) 


( 2 - 1 - 2 ) 


where f is the input specific force and f is the sensed specific force. Xl, Xs, and Xa are 
time constants and K is also a constant. Xl = 5.33 sec, Xs = 0.6 sec. Xa= 13.2 sec, K = 0.4. 
2.2.Reference Frames 

In describing the development of motion cueing algorithms it is convenient to 
employ several reference frames. These reference frames are defined below and are 
shown in Figure 2. 1 . 
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2.2.1. Frame Fr$ 

The simulator reference frame Frs has its origin at the centroid of the simulator 
payload platform, i.e. the centroid of the simulator's upper bearing attachment points. It 
is fixed with respect to the simulator payload platform. Xs points forward and Zs 
downward with respect to the simulator cockpit. Ys points toward the pilot's right hand 
side. The x-y plane is parallel to the floor of the cockpit. 

2.2.2. Frame FrA 

The aircraft reference fimne FrA has its origin at the same relative cockpit location 
as the simulator reference fimne Frs. FrA has the same orientation for Xa, Ya, and Za 
with respect to the cocl^it as the simulator frame Frg. 

2.2.3. Frame Free 

The aircraft center of gravity reference frame Free has its origin at the center of 
gravity of the aircraft. Frame Free has an orientation for Xcg, Ycg, and Zcg that is 
parallel to reference fiames Frs and FrA. 

2.2.4. Frame Frps 

This is a frame attached to the simulator pilot's head with its origin located at a 
point midway between his left and right vestibular systems. Frame Frps hag an 
orientation for Xp$, Yps, and Zp$ which is parallel to FrA and Fr$ . 

2.2.5. Frame FrpA 

This is a frame attached to the aircraft pilot's head with its origm located at a point 
midway between his left and right vestibular systems. Frame FrpA has the same 
orientation for Xpa, Ypa, and Zpa as Frps . 
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2.2.6. Frame Fri 

The inertial reference frame Fri is earth-fixed with Zi aligned with the gravity 
vector g . Its origin is located at the center of the fixed platform motion base. Xi points 
forward and Yi points to the right hand side with respect to the simulator pilot. 

2.2.7. Reference Frame Locations 

In Figure 2.1 are four vectors which define the relative location of the reference 
fi:ames. Ri defines the location of Frs with respect to Frj. Bs defines the location of Frps 
with respect to Frs. Similarly, Ea defines the location of FrpA with respect to FrA- Ecg 
defines the location of FrA with respect to Free- 
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2.3. Coordinate Transformation 

Consider a body having both translational and rotational motion with respect to 
Fri. In Figtire 2.2, the inertial translational acceleration of a point b on the body located 

at a distance 5:=8,i + 8J + 5,k from the origin of any body reference frame Fre can be 
expressed by: 

a[ =aB+6b+2©Bx5®+®Bx8®+coBx(©^x6b) (2-3-1) 

v^iiere 

§3 is the acceleration of the origin of Fre with respect to Frj. 

©B is the angular rate of Frs with respect to Fn. 

p,q,r are the three components of © ^ , i.e., ©B=pi+qj + rk. 

** B * B 

Assuming point b is fixed with respect to Fre, then 8^=6^= 0. 

al = 4 +©'b x8® +©L x(©L x8®) (2-3-2) 

= aL + [-5,(q^ + r^) + 8,(pq-f)+8,(pr + ^] i 
+ [ 5Jpq + f) -8y(p^ + r^) + 8^(qr -p)]j 


+ [ 5,(pr-q) +8y(qr + p)-8,(p^+q^)]k 



Figure 2,2. Coordinate System for Inertial Acceleration. 
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2.3.1. Euler Angles 


The orientation between two reference frames Fri and Fre can be specified by 


three Euler angles: p = [(j» 0 which define a sequence of rotations which will cany 
Fri into Fre . A vector Y expressed in the two frames can be related by the transformation 


matrix L® ( Frj -> Fre ) or Lbi ( Fre — > Frj ): 
Y® = LbY’ and Y' = LbiY® 

where 


Lib Lgj L^j 


(2-3-3) 


Lbi — 


cosGcosy 

cosBsiny 

-sin0 


sin (|> sin 0 cosy -cos <j) sin v|/ cos<j»sin0cosy -i-sin(|)siny 
sin<}>sin0siny + cos<j>cosy cos <j> sin 0 sin y — sin <1> cosy 
sin(j>cos0 cos<|>cos0 


The angular velocity of Fra with respect to Fri can be related to the Euler angle 


rates P by the following. Let represent the components of this angular velocity in 
frame Fre, then 


P — Tg o B 


where 


and 

where 



1 sin<j>tan0 cos<|>tan0 

Tb = 

0 cos(|> -sin<|) 


0 sin«|)sec0 cos^secG 


Tb‘P 


1 0 -sin0 

Tb* = 

0 cos(|> sin<|>cos0 


0 — sin«|) cos(|)cos0 


(2-3-4) 
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F^re 23. Illustration of a Six Post Synergistic Motion Base Geometry. 
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2.4. Actuator Extensions 


The geometry of a six degree-of-ffeedom synergistic motion base is given in 
Figure 2.3. The relevant vectors relating the locations of the upper emd lower bearings of 
the j-th actuator are shown in Figure 2.4: 



Figure 2.4. Vectors for the j-th actuator. 

In Figure 2.4 Os and Oi are the centers of the payload platform and fixed platform 
respectively. Os and Oi are respectively the origins for Frs and Fri. 

It can be seen that the relation among those vectors is: 


r,+a]=Rj=b;+«j 

Then the actuator length vector can be foimd from: 

Expression of f j in Fri is desired: 

iS=A;+R,-B; 


(2-4-1) 


(2-4-2) 


(2-4-3) 
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Where Aj are the coordinates of the upper bearing attachment point of the j-th actuator 

in Fr* and B] are the coordinates of the lower bearing attachment point of the j-th actuator 
in Fri. 

The actuator extensions can then be found from 


= (Ls,(t)-Ls,(0))Af +(S'(t)-S'(0)) 


(2-4-4) 


Usually Ls,(0) = I and S'(0) = 0, where I is the identity matrix. Therefore 


= ALs,Aj + AS’ 


(2-4-5) 


Using small angle approximation AL^, can be expressed as: 




0 ^ G. 


Vs 0 


~^s 4*S 0 


(2-4-6) 


Employing Equations (2-4-5) and (2-4-6) the actuator extensions can then be calculated. 

It is observed that a smaller A] will result in smaller actuator extensions for a given 
simulator rotation angle. This mformation might be useful for simulator hardware 


2.5. Input Scaling and Limiting 

Scaling and limiting are applied to both aircraft translational input signals a^ and 

rotational input signals ©^. Scaling and limiting modify the amplitude of input 
uniformly across all frequencies. Limiting is a nonlinear process that clips the signal so 
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that it is limited to be less than a preselected magnitude. Scaling and limiting can be used 
to reduce the motion response of a flight simulator. Two input scaling and limiting 
algorithms were used in the current simulation software. They were suggested and used 
in [3], [4], and [5]. 

2.5.1. Linear Input Scaling in Combination with an Input Limiting 

The first algorithm is characterized by a linear input scaling in combination with 

an input limiting. Each component of a^ and in different degrees of fi^edom is 
scaled and limited separately but in the same manner. The scaling and limiting of aAx, the 
X component of a^ is given as an example, where Sx is the slope from -Xi to Xj: 


^x^Ax 

I^Axol— X] 

s.x, 

^Ax - Xj 

-SxX, 

^Ax^-X, 


(2-5-1) 



INPUT 


2.5.2. Nonlinear Input Scaling 

The second algorithm is characterized by a nonlinear input scaling. Input limiting 
is not used. Each component of a^ and in different degrees of freedom is scaled 
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sqjarately but in the same manner. The scaling of bax, which is the x component of a 
is used as an example: 


a 


X 


Sx^xo laxol^X, 

SxSxo - 0.7S,(a,o - X, ) a,o > X, (2-5-2) 

Sxaxo-0.7S,(a,o + X,) a,„<-X, 



Figure 2.5.2, Nonlinear Input Scaling. 
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3. Structure of Simulation Software 


The execution of the software begins with input of a set of simulation commands. 
The simulation commands contain selection of input degree of ft^eedom, input type, input 
magnitude and duration, duration of simulation, and type of cueing algorithm. Each of 
the input channels for six degrees of fieedom can be selected. Maximum duration of 
simulation is set to be 40 seconds. There are six aircraft input options: step, pulse, pulse 
doublet, sinusoidal, half sine pulse, and ramp to step. Aircraft translational acceleration 

at the point corresponding to the centroid of the simulator payload platform a^^ and 

aircraft angular acceleration © ^ are used as input vectors. 

From the aircraft inputs the aircraft response and the simulator response are both 
calculated. Next the aircraft pilot's sensation and the simulator pilot's sensation are 
calculated and compared. Actuator extensions are generated based on the simulator 
responses output by the cueing algorithm. 

The aircraft response is assumed to follow the input conunand without error. 

Then the acceleration at the pilot's head §p^ can be calculated from a^ and By 
subtracting the gravity vector g the specific force on the pilot's head f ^ can next be 

calculated. By passing and f^ through the vestibular model the aircraft pilot’s 
sensation will be generated. 

The simulator response is calculated by passing the aircraft commands through 
the selected cueing algorithm and the platform dynamics filter. The cueing algorithms 
are the kernels of this software. They are responsible for maximizing the motion cueing 
effects while restricting the physical motion to be within the displacement, velocity and 
acceleration capacity of the motion system hardware. The cueing algorithm outputs the 
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desired translational and rotational platform positions S^and P^wdiich are used to 

compose the desired actuator commands f ^ that will drive the simulator platform. 

In a real simulation, the platform dynamics will cause some error between the 
platform motion and the motion commands. Passing the motion commands through a 
filter that is a model of the platform dynamics can simulate this. Then the filter outputs 

P and $' are assumed to be the real platform positions. Based on p and S' both the 

— s 

simulator translational acceleration a^ and angular rate may be obtained. The 

specific force on the simulator pilot's head f ^ is also available. By passing ©j and f ^ 
through the vestibular model both the sensed angular rate and sensed specific force are 
obtained. The flowchart of the simulation software is shown in Figure 3.1. 
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SPECinC FORCE SENSED 
BY SIMULATOR PILOT ) 


Figure 3.1. Overall Structure of the Software 
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4. Descriptions of the Four Washout Algorithms 

Four algorithms are used in this study. The first is known as the classical 
algorithm. This type of algorithm is generally denoted in the literature as a linear cueing 
algorithm [3], [4]. The second algorithm evaluated in this study is the NASA Langley 
Research Center developed "Coordinated Adaptive Washout Algorithm" [5]. The next 
cueing algorithm reviewed in the study is a variation of the Coordinated Adaptive 
Washout Algorithm. This was developed by L. Reid and M. Nahon at the University of 
Toronto [3], [4]. The Optimal algorithm, the fourth algorithm employed in this study is 
that developed at MIT by Sivan, et al. [6] and implemented as described in [3] and [4]. 

The basic task of the washout algorithms is to create a specific force vector and an 
angular velocity vector at the pilot's location in the simulator approximating those that the 
pilot would experience in an actual aircrafi. The translational and rotational motion 
effects on the simulator pilot are expected to approximate the motion effects on the 
aircrafi pilot: 

■pS ^ x*A S ^ 

fps ~ i-PA ’ ^PS ~ ^PA 

The relation between the specific force acting on the simulator (aircrafi) pilot and the 
specific force at the origin of the simulator (aircrafi) frame can be found fi^om Equation 
(2-3-2); 

fS „s „s 
£ps ~ ^ps 8 

= a® + ® s X Rj + X (®s X Rs ) - g^ 

= fs +“s ^ Es Es) 

fpA=a^A-g^ 

= ® A ^ Ea + ® A A X Ea ) - s'" 

= f A + ® A ^ Ea + ® A (® A ^ Ea ) 
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Usually Ra ® s “ “ 5 ^ra • Thus the washout algorithms attempt to 

achieve: 


f s ^ f A „A „f 

fs^fA =aA"g 


®S*®A 


f CO* or a^ , CO* are used as washout algorithm inputs . 

4.1. Classical Algorithm 

This algorithm employs aircraft body axes acceleration a* and angular velocity 


CO* as the aircraft state vector elements that provide the input to the cueing algorithm. A 
linear scaling in combination with a limiting as described in Section 2.5 is applied within 
the algorithm to modify the input The architecture of the classical approach is such that 
there are separate filters for the translational degrees of fiuedom and the rotational 
degrees of fi-eedom with a crossover path to provide the steady state or gravity align cues. 
This algorithm behaves like open-loop control. The details of this algorithm are 
presented below. 

4.1.1. Translational Degrees of Freedom 

The aircraft acceleration vector a^ is first scaled and limited. This scaling can be 
either linear or nonlinear. It should be noted at this point that it is not the scaling that 
makes the cueing algorithm either linear or nonlinear, but rather it is the formulation of 
the washout filters that is responsible for that char^teristic. For the classical algorithm a 
scaling and limiting scheme as described in Section 2.5.1 is used. After scaling and 

l imit ing, the aircraft specific force f^ is computed and then transformed fixrm the 
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simulator frame Fr$ into the inertial frame Fi^, whereupon the inertial frame vector 


subtracts the gravity vector, with the difference f j filtered by a high pass filter (transfer 
function) of the form 

s^ 

TF = (4-1-1) 

(s^+2q©„s + (0^)(s+<Db) 

Where On is a second order system imdamped natural frequency and o>b is a first order 


system break frequency, q is a second order damping ratio. The output is then 
integrated twice to provide the simulator translational position S*- 
4. 1 .2. Rotational Degrees of Freedom 

The rotational rate vector is first scaled and limited by the same scheme for the 
translational channel. The resulting vector is transfonned to the Euler angular rate. The 
Euler angular rate is then filtered through a high-pass filter (transfer function) of the form 


TF = 


(S + fflJ^ 


(4-1-2) 


The output is then integrated to provide the desired simulator Euler angular position 
corresponding to the aircraft rotational input. 

4.1.3. Tilt Coordination 

Sustained translational acceleration is sensed by the pilot as a long term change in 
the magnitude and direction of the specific force in the absence of rotational motion. 
This cannot in general be simulated by translatioiuil motion due to motion-base travel 
limits. It is possible to alter the direction of the steady-state specific force experienced by 
the pilot in the simulator by tilting the cab. It has become common practice in flight 
simulators to employ cab tilt to simulate the effect of sustained translational inertial 


31 


acceleration. Because this tilt coordination process cannot alter the long term magnitude 
of the specific force vector, it is an approximation to the desired effect. [3], [4] 

The formulation of the tilt coordination Euler angles starts with a^ as shown 


in Figure 4-1-1 . a^ is formed by filtering a* through a low pass filter of the form: 


TF = 


0 ). 


s^ +2qto„s + (o^ 


(4-1-3) 


Cl can be considered to be the specific force components that are to be simulated 
through simulator cab tilt. In the absence of other cab rotational displacement and 
translational motion, f^ in FrA is shown in Figure 4.1.1 (a). By tilting the simulator cab 

f| can be rotated with respect to Frs (Figure 4.1.1 (b)). It is found that this change in 
the direction of the sensed sustained specific force by the simulator pilot is quite useful in 
the simulation of low fiequency aircraft acceleration components. When the desired tilt 

angle a , = a o , f® = k f ^ , where k = g / 1 f ^ | , f| is successfiilly aligned. It is easy to 
see that the desired tilt angles a , are: 


‘I'sL = tan"' (-fry > g) « > g 

esL = tan''(f,^/g)«f,5c/g 
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F^;ure 4.1.1. Simulator Cab Tilt 

P is usually restricted under the threshold of pilot perception for avoiding undesired 

—ST 

cue. 

4.1.4. Summation of Two Rotational Channels 

The summation of P and P will yield P , the angular position of the 
simulator. Lsi and Ts can be formed by Equations (2-3-3) and (2-3-4). Then the 
simulator translational position S‘ and the angular position p can be transformed from 

degree-of-freedom space to actuator space. These are the actuator lengths required to 
achieve the desired platform translation. 

The block diagram for the classical algorithm is shown in Figure (4-1-2). 


33 




Figure 4.1.2. Classical Algorithm. 









4.2. NASA Coordinated Adaptive Washout Algorithm 

This algorithm was developed at the NASA Langley Research Center in 1977. 
The intent of the algorithm is to adapt the severity of the simulator washout filters 
according to the current state of the simulator. In this way it should be possible to make 
full use of the simulator motion system at all time. This algorithm also employs 

a* and as the input to the cueing algorithm. A nonlinear scaling as described in 
Section 2.5 is applied to modify the input. The general architecture is similar to the 
classical algorithm in the aspect that there are separate filtering channels for the 
translational degrees of freedom and the rotational degrees of fiwdom with a crossover 
path to provide the steady state or gravity align cues. The block diagram for this type of 

adaptive algorithm is shown in Figure 4-2. The signals fj and represent the inertial 

frame acceleration components and Euler angles which if t^plied to the simulator frame 
Frs will produce specific force and rotational velocity components in Fr$ identical to the 

corresponding FrA components in the aircraft. P^ is passed through the rotational channel 
with an adaptive gain 6 to produce a simulator angular rate command. The gravity vector 
g is added to f j to yield a new specific force f\ that will actually be simulated as 

simulator motion, f ^ is passed through a translational channel with an adaptive gain X, to 
produce a simulator translational acceleration command and passed through a crossover 
tilt channel with a fixed gain y to generate a simulator angular rate command. These 
commands drive the simulator to the desired translational and angular positions. 

S* and S* are employed as feedback as shown in Figure 4-2. Adaptive parameters 
X and 5 are continuously adjusted according to the current state of the simulator and 
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aircraft input. The control equations for this cueing algorithm are presented below. The 
motion equations are separated and dealt with as four parallel modes: pitch/surge, 
roll/sway, yaw, and heave. 

4.2.1. Pitch/Surge Mode 
Control law: 

“ ^xfx ~^X^x ~®X^x (^*2-1) 

®s ~ Yxfx "*'^X®A 

whCTe dx,Cx & Yx are fixed parameters for the pitch/surge mode; anddx are the 
pitch/surge adaptive parameters which are continually adjusted in an attempt to minimize 
the instantaneous value of the cost function. 

Cost ftmction: 


h = -X)" +^(e* -9,)’ (4-2-2) 

where Wx, bx, and Cx are constant weights which penalize the difference in response 
between the aircraft and simulator, as well as restraining the translational velocity and 
displacement in the simulator. 

Steepest descent for the adaptive parameters: 






(4-2-3) 


where , and K;g^ are constants. The first right hand side term of each 

equation in (4-2-3) defines that the change of the adaptive parameter (Xx or 5*) is toward 
a minimum cost function and also defines the rate of change together with the second 
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right hand side term. The second right hand side term also restrains the deviation of 
either Xx or 5x from their original values. 

4.2.2. Roll/Sway Mode 

The control law, cost function and steepest descent for the ad^tive parameters 
are in the same form as for the pitch/surge case. Therefore by a substitution of y for x 
and <|> for 0 the roll/sway motion equations may be obtained. Note that the adaptive 
parameters, cost function weights and steepest descent constants are imique for the 
roll/sway mode with a subscript y replacing the subscript x. 

4.2.3. Yaw Mode 
Control law: 

¥ = (4-2-4) 

where e^ is a fixed parameter while is an adaptive parameter. 

Cost function; 

K =|(Va -Vs)^+^Vs (4-2-5) 

viiere b^, is a constant weight in the cost function. 
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Steepest descent: 


where K„ and K=„ are constants. 

My “My 

4.2.4. Heave Mode 
Control law: 

Sl=Ti,f^-d,SL-e,S^ 

where y\^ is a fixed parameter while dz and ez are adaptive parameters. 
Cost function: 

1 , , b, , C, 

Jz = 2^fl-z)'+YZ'+Y^ 

where bz and Cz are constant weights in the cost function. 

Steepest descent: 

|~ + Ki^(nzo - Tlz) 

where K„ and K;„ are constants. 

»lz 


(4-2-6) 


(4-2-7) 


(4-2-8) 


(4-2-9) 
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4.3. UTIAS Coordinated Adaptive Algorithm 

This algorithm was developed at the University of Toronto in 1985. The general 
philosophy and structure of this algorithm are the same as those of the NASA adaptive 
algorithm. The UTIAS adaptive algorithm adapts the severity of the simulator washout 

filters according to the current state of the simulator, a^ and ® ^ are the input to the 
cueing algorithm. A linear scaling in combination with a limiting as described in Section 
2.5 is applied to modify the input. The block diagram is shown in Figme 4-3. The 

signals f} and represent the inertial frame acceleration components and Euler angle 

rate \^ch if applied to the simulator fiame Frs will produce specific force and rotational 
velocity components in Frs identical to the corresponding FrA components in the aircraft. 

is passed through the rotational channel with adaptive gain p *3 to produce a simulator 

angular rate command, f^ is passed through a translational channel vrith an adaptive 
gain pxi to produce a simulator translational acceleration command and passed through a 
crossover tilt channel with an adaptive gain px 2 to generate a simulator angular rate 
command. These commands drive the simulator to the desired translational and angular 

positions. S* and S* are used as feedback as shown in Figure 4-3. Ad^tive parameters 
Pxi, Px2 and px 3 are continuously adjusted according to the current state of the simulator 
and aircraft input. The difference between the UTIAS and the NASA adaptive 
algorithms is that the UTIAS algorithm made more pa.ameters adjustable and employed 
more terms in the cost function in an attempt to improve the effect of the ad^tive 
washout filter. The parameter in the crossover tilt channel is adjustable in the UTIAS 

algorithm and remains fixed in the NASA algorithm. The absolute values of p and P 


40 



along with the deviation of the adaptive parameters from their initial values are employed 
in the cost fimction. The control equations for this cueing algorithm are presented below. 
The motion equations are separated and dealt with as four parallel modes: pitch/sxirge, 
roll/sway, yaw, and heave. 

4.3,1. Pitch/Surge Mode 
Control law: 

Sx = PxifJx-kxiS‘-k:0!Si (4-3-1) 

Oj = lim(px2 fix) + Px3 

where kxi and kx 2 are fixed parameters and parameters pxi, px 2 , and Px3 are adaptive 
parameters which are continually adjusted in an attempt to minimize the instantaneous 
value of the cost fimction. 

Cost fimction: 

Jx= 0.5 [ y, (f|x-Si)2 + Wx,(0^-0s)" 

+ Px (Wx2 + Wx3 S',' + Wx4 0^ + Wx5 0s) (4-3-2) 

+ Wx6 (Pxl - Pxio)^ + Wx7 (Px2 - Px20)^ + W*8 (Px3 " Px3o)^ ] 

where yx , Px , and Wxi are constant weights. In (4-3-2) the first group of terms in the 
cost fimction penalizes the difference between the aircraft and simulator responses. The 
second group restrains the translational velocity and displacement of the simulator. The 
third group restrains the deviation of the adaptive parameters fiom their original values. 
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4.3.2. Roll/Sway Mode 


The control law and cost function are in the same form as for the pitch/surge case. 
Therefore by substituting y for x and ^ for 0 the roll/sway equations may be obtained. 

4.3.3. Yaw Mode 
Control law: 

Vs = Py Va - jvsdt - K 2 Vs (^-3-3) 

where k^, and parameters and p^ is an adaptive parameter. 

Cost function: 

Jy = 0-5 [(Va- Vs)^+Py (Wm-I Vs^ + Wy2\|/s) + Wys(py - pyo)'] 

(4-3-4) 


where and Wyi are constant penalty weights. 

4.3.4. Heave Mode 
Control law: 

s'. = P. £lz - k., Js'.dt - S'. - k„ s; (4-3-5) 

where k^ are fixed parameters and pz is an adaptive parameter. 

Cost function: 


h = 0-5 [(flz - Kf + Pz(Wz,S-^ + + W^(p,-p^)^] 


(4-3-6) 


where pz and Wzi are constant penalty weights. 
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4.3.5. Steepest Descent for the Adaptive Parameters 
For all four channels the steepest descent has the same form: 





" 5p, 


(4-3-7) 


With the subscript c corresponding to x, y, and z, for \|/ the subscript] corresponds to 1, 
2, and 3, and Gc are constants. 
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Figure 4.3. UTIAS Adaptive Washout Algorithm. 













4.4. MIT/UTIAS Optimal Washout Algorithm 
4.4.1. Optimal Control Theory 

A linear optimal control theory can be applied to the aircraft simulation problem. 
The definitions and solutions of linear optimal control problems have some simple 
standard forms, which make it convenient to apply the optimal control theory and find the 
solution to a specific optimal control problem. 

I. Deterministic Linear Optimal Regulator Problem 

The problem can be illustrated by Figure 4.4.1. The problem is to determine 
some constant matrix F that relates the control system input u to the system states x and 
minimize s the optimal criterion J. 



Figure 4-4-1. Deterministic Linear Optimal Regulator. 
System equations: 


fx = Ax + Bu 


< ^ with X = x(0) at t = 0 . 

[ y = Cx 

(4-4-1) 

Criterion: J= x + u’^Rj u]dt + x^(t,)P, x[t,) 

(4-4-2) 


The problem of determining an input n for which J is minimized is called the 
deterministic linear optimal regulator problem. 
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II. Stochastic Linear Optimal Regulator and Tracking Problems 

(a) Regulator Problems with Disturbances: the Stochastic Regulator Problem 

The problem can be illustrated by Figure 4.4.2. Disturbance input y is a filtered 
white noise with limited bandwidth. The problem is to determine a constant matrix F that 
determines u by relating it to the system states x and the disturbance input y so that the 
optimal criterion J is minimized. 


w 

V 



FI, F2 are partitions 
of F. 

F = [F1 F2] 


Figure 4.4Jt. Stochastic Linear Optimal Regulator. 


System equations: 


r X = Ax + Bu + V 



y = Cx 

= AdisHdi, + w 


V = C^isXdis 


is white noise. 


(4-4-3) 


Define x = 


X 




, then 
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(4-4-4) 

(4-4-5) 


~ P 0 

where Ri = Ra C , Pj = J ^ , Ri, R 2 , R 3 , and Pi are constant matrices, and E is 
defined as the mathematical mean of a statistical variable. 


(b) Stochastic tracking problems. 

This problem can be illustrated by Figure 4-4-3. 2Jr is the state of a reference 
system which has a white noise input, s is the state of the control system which strives 
for following the reference system. The problem is to determine a constant matrix F 
which determines u by relating it to the state of the reference system ^ and the state of 
the control system x so that the optimal criterion J is minimized. 

FI, F2 are partitions 
of F. 

F = [F1 F2] 

Figure 4.43. Stochastic Tracking Problem 
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System equations: 


X = Ax + Bu 
y = Cx 

< . “ . w is white noise. 

Xr = A^Xr +5^ 


Define x = 


X 


, then 




A 0 ' 


B‘ 


o' 


x-h 


U + 


.0 


0 


w 


y=y-y,=[c -c,]x 

Optima] criterion: 

j = E I llcy - Rj(y - y,) + m] dt | 

= e| |*'[x’^R,x + u’^R2uj dt 
where Ri=[C -Cr]^R 3 [C -CrJ 


(4-4-6) 


(4-4-7) 


(4-4-8) 


Equations (4-4-6) and (4-4-8) can be generalized in a common standard form: 


System equation: 


[^x = Ax + Bu + Hw 
1 y = Cx 


Optimal criterion: J = E | [x’^R, x + u’^Rj u]dt + x^P, x | 


(4-4-9) 


(4-4-10) 


where Ri and R 2 are constant matrices and E is defined as the mathematical mean of a 
statistical variable. 
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in. Solution to the Optimal Control Problems 

There is a common linear optimal solution for the deterministic regulator 
problem, the stochastic regulator problem and the stochastic tracking problem. The 
solution has a standard form: 


U = -Fx (4-4-11) 

where 

F=R-*B^P(t) (4-4-12) 

and P(t) is the solution to the algebraic Riccati equation: 

-P = R, -PBRjB^P PA + A^P (4-4-13) 

with boundary conditions x = x(0) at t = 0 and P(ti) = Pi. 

When ti approaches infinity, it is proven that P has a steady-state solution which satisfies 
the equation: 

0 = R,-PBRjB’^P + PA -h A^P (4-4-14) 

4.4.2. Aircraft Simulation Problem Definition 

The problem is to determine a linear filter matrix W(s) wdiich relates the simulator 
motion states us to the aircraft motion states Ua so that some cost function or criterion 
which constrains the pilot’s sensation error e and the simulator motion simultaneously 
will be minimized. The structure of the problem is illustrated in Figure 4.4.4. 
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Figure 4.4.4. Aircraft Simulation Problem Structure 
Since the control strives for tracking the output of the aircraft pilot’s vestibular 
system, the problem is most likely to be treated as a tracking problem. The aircraft 
motion can be quite variable. Therefore, it is reasonable to use a filtered white noise that 
contains sufficient frequency components to represent the aircraft motion states. 
Therefore, the problem is a stochastic tracking problem. 

4.4.3. MITAJTIAS Development of Washout Filter Coefficients 

The optimal algorithm documented in this section was developed at the 
University of Toronto in 1985. The optimal washout filter design problem is formulated 
as follows. The actual aircraft and simulator sensation systems SensationA and 
Sensations are given along with the input signal Ua that drives SensationA. A properly 
constrained operator W(s) can then be found that generates the simulator input Us to 
Sensations on the basis of the input to SensationA, such that the pilot senation error £ is as 
small as possible. A mathematical model of the vestibular organ is used. The 
optimization criterion that is selected is the mezm-square difference between the 
physiological outputs of the vestibular organs for the pilot in the aircraft and for the pilot 
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in the simulator. The structure for the optimal filter is assumed given and the task is to 
find the optimal values for the parameters of the filter. 

The optimal algorithm in this section generates the optimized transfer function 
W(s) by an off-line program. W(s) is then applied on-line to generate simulator motion 
commands. 

W(s) relates the aircraft pilot sensation input to the simulator pilot sensation input 
that is assumed to be approximately identical to the simulator motion base input 

Us - W(s) ua , where ua consists of aircraft body axes acceleration and angular 


displacement . us will be used as the command to drive the motion base. A linear 

scaling in combination with a limiting as described in Section 2.5 is applied to modify the 
input. 

W(s) is optimized by minimizing the cost function a = J J dt , where 
J = E|e"Qe + p(ujRUs + ^R.yJ} 


£ = Xs-Xa 

where ^ and y, are the pilot's sensations in the aircraft and simulator environments. 

For uncoupled system equations, yjJ has different meanings for different degrees of 
fieedom: 


For pitch/surge equations: 

- [ JJu>= 

Jusdt “si 

For roll/sway equations: 


Jvsdt v‘] 

For yaw equations: 

&?-[ JJv 

Jvsdt] 

For heave equations: 

- [ JK*' 

Jw'dt w^ 
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Q and are positive semi-definite matrices, R is a positive definite matrix and p is a 
positive scalar. These parameters determine the weight of each term in the cost fimction. 

The optimal washout is performed in the pilot head frames FrpA and Frps. It is 
claimed in the UTIAS report [3], [4] that this frame selection will avoid sensation cross- 
coupling where the sensation cross-coupling may make the calculation of W(s) more 
complicated. 

The washout filter coefficient W(s) will now be generated. The sensation system 
for the pilot in both the aircraft and simulator environments are given as: 


Sensation^: 

Xa = A, Xa+Bi Ha 
y^=c, Xa+D, Ua 

(4-4-15) 

Sensations: 

Xs = A, Xs+Bi Hs 

ys=C, Xs+D, Us 

(4-4-16) 


Where it has been assumed that the same motion sensing system dynamics can be applied 
both in the aircraft and in the simulator and all system matrices are taken to be time- 
invariant. 

Assume consists of filtered white noise: 


x„=A. x„+B. w 


Ha =x„ 


(4-4-17) 


Adjoining all the above systems, the optimal controller equations may be generated: 


X = AX + Bus+Hw 
Y = CX + Dus 


(4-4-18) 


with the cost function a 




Jdt, 


\vdiere J = eIy’^G Y + p UgRus} 


(4-4-19) 
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o' 
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nl 
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B, 


UiJ 
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. 0 » 


0 


c, c, 


Q c. 


-D. 

0 


D = 


D, 

0 


Q Q 

0 Rd 


H = 


0 

0 


where p and R are positive scalars. 


By matrix substitution and manipulation Equations (4-4-18) and (4-4-19) transform to the 


standard form of the stochastic linear optimal regulator problem: 


X = A'x + B u! + H w 

J = E I £ [x^R;x + if'^Riu:] dt I 

vvdiose solution was given in Section 4.4.1 : 

u! = -F, x> 

where 

F, = r;‘b^p 

and P is the solution of the matrix Riccati equation 


(4-4-20) 


(4-4-21) 


0 = R; -PBRj'B^P + A'^P + PA' 


(4-4-22) 
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with definitions of the new notations used in the above equations; 


R, = G C ; Ri2 = G D ; R 2 = R + G D; 


R; = R, - RijR^'R^j ; u! = Us + R 2 RJ 2 X ; (4-4-23) 

A' = A-BR'‘R[j ;x = X 
Now Us and x can be related: 

Us = u!-R'‘R^ 2X = -R-' (B^P + R^j) x= -Fx (4-4-24) 


where -F = -R;’ (B^P + RI) . 
Partition -F into - [FI F2 F3], then 


Us = - [ FI F2 F3 ] 



Taking the LaPlace transform on (4.4.1) and (4.4.3) yields: 


(4-4-25) 


X.=(sI-A)-'Bu^ 

Xs=(sI-A)-'BUs 


(4-4-26) 


By the substitution of (4-4-26) into the LaPlace transform of (4.4.25), the relation 
between Us(s) and iia(s) is finally foimd: 


Us(s) = W(s) u,(s) 

where (4-4-27) 

W(s) = [-I + F2(sI-A + BF2) *B] (FI (si- A)*B + F3[ 

W(s) is the optimized open-loop transfer function linking us(s) and Ua(s), the optimal 

algorithm controller implemented in the UTIAS report. The block diagram for the 

optimal algorithm is shown in Figure 4.4.5. 
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LIMITING 



Figure 4.4.5. Optimal Washout Algorithm. 













This Page Intentionally Left Blank. 


56 



5. Changes Made in Project to Original Algorithms 

For evaluation purposes, changes were made to the four algorithms. The NASA 
adaptive algorithm and the three washout algorithms implemented by UTIAS (classical, 
adaptive, and optimal) are accommodated in one program as described in Section 3. 
Some of the subroutines in the NASA software have functions overlapping with the 
functions of some subroutines in the UTIAS software. These overlapping subroutines are 
not included in the current evaluation software. 

Limits on angular tilt rates were not included in the optimal and the NASA 
adaptive algorithm. In the current software these limits are included. The parameters for 
scaling translational and rotational inputs were different for each algorithm. For this 
project a scale factor of 1.0 is used in all algorithms for convenience of comparison 
between the aircraft motion and the simulator response. 

When either the UTIAS adaptive or NASA adaptive algorithm was run, 
convergence to steady state oscillations was observed under rotational input. The 
adaptive algorithms were modified in such a manner as to influence the response of the 
simulator corresponding to the aircraft rotational input. The modification does not affect 
the response of the simulator to the aircraft translational input. The difference is that the 
simulator will perform pure rotation to simulate an aircraft pure rotational input. The 
original algorithm generated some translational acceleration under aircraft pure rotational 
input and this translational acceleration is occasionally unstable. The rotational response 
of the simulator to rotational input was not affected. This topic is further discussed in 
Section 5.1. 
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When either the UTIAS adaptive or NASA adaptive algorithm was run, unwanted 
spikes on the specific force of the pilot's head occurred whenever the input contains some 
sharp change of translational acceleration. These spikes were in the opposite direction to 
the correct response direction. The UTIAS report [3]; [4] mentioned these unexpected 
spikes but did not explain the correct reason for their occurrence. No effort was made to 
eliminate or decrease the spikes. The cause for the spikes was identified to be angular 
acceleration when the simulator is tilting for compensating aircraft translational 
acceleration input. This angular acceleration can be so large that it can overpower all 
other simulator motions at some points in time and give the pilot a perception in the 
opposite direction to which is expected. The angular acceleration is limited in the current 
software. The spikes were effectively reduced and the overall response of the simulator 
did not change significantly. This topic is further discussed in Section 5.2. 

A platform dynamics filter was added. A model of the actual platform dynamics 

©n* 

is not currently available. A transfer flmction "2 is used to q>proximate 

S + 2^©n + ffin 

the platform dynamics, where q = 0.707, ©n = 10 7t rad/s. This filter represents the 
specifications for the real platform design. In the output plots in Appendix A, the 
notation 'desired simulator displacement' and 'desired simulator angular position' means 
the ou^ut of a washout filter. This ou^ut has not yet passed through the platform 
dynamics filter. The notation 'actual simulator displacement' and 'actual simulator 
angular position' means the output that has passed through the platform dynamics filter. 
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5. 1 . Instability in the UTIAS Adaptive Algorithm 

In some rotational input cases, the UTIAS adaptive algorithm was unstable. In 
the UTIAS report, Volume 2, [4] the instability of the adaptive algorithm was mentioned 
but the reason was not discussed. The report suggested restricting the adaptive 
parameters to eliminate the instability. But this restriction will make the algorithm lose 
some of its adaptive characteristics, i.e., the adaptive algorithm becomes less 'adaptive', 
and this restriction may not eliminate the instability completely in some cases. It is 
necessary to find the cause of the instability before trying to eliminate it. 

As shown in Figure 4.3, the algorithm tries to duplicate f^ in the simulator 
frame. The simulator rotates at an angle 6s to simulate an aircrafi rotation of 6 a- It is 

always true that |0 s|<|9a| unless both are zero. In the pure rotational input case in 
which §A = 0, the aircraft and simulator angular positions are plotted in Figure 5. 1 below; 



Figure 5.1. Aircraft and Simulator Angular Position under Pure Rotational Input 

It is observed that with the doublet angular acceleration input, both the aircraft 
and the simulator angular positions reach steady state after some time. In these steady 
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states, g imposes a force on the pilot's head in the aircraft and a force g® on the 
pilot's head in the simulator, where g"' g® . For example, in the pitch test case, 

g'^= [-g* sin 0^ 0 g^cosG^f 


g®= [ -g * sin 9g 0 g * cos 0s 

with f^ = g^- g® = [-g*(sin 0s -sin 0^ ) 0 g * (cos Gg - cos 0^) ]^ . 


(5-1) 


The UTIAS adaptive algorithm attempted to compensate for f by passing S* and P 

—ST 

through both the translational and crossover tilt channels. Now look at the formula for 
generating the simulator translational acceleration; 

S' =Pxi*fL-kx.*S^ k.2* S* (5-2) 


where kxi and kx 2 are constant scalars and pxi is an adaptive parameter, f^ is not zero in 

steady state since |Gs| < |Ga| Since all washout algorithms should attempt to wash out 

s' to zero in steady states, the adaptive algorithms would attempt to bring the simulator 
to its neutral position in the steady state. When the simulator reaches its neutral position, 

i.e., s' = 0, it holds that S^ + kx 2 * S' = pxi *' f ^ - Since f^ is not equal to zero, S* and 

S^ cannot both be zero. Then the simulator cannot stay in the neutral position, i.e. it must 
continue moving. Because the adaptive algorithms are attempting to return the simulator 
to its neutral position in the steady state and the simulator cannot stay in its n^tral 
position, the simulator oscillates around its neutral position. This is the cause of the 
instability under rotational input This analysis is consistent with the results generated by 
the original adaptive algorithms before revisions were made. 
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The difference between g'^and g® under rotational input in fact cannot be 
eliminated. It cannot be eliminated by simulator rotation since 1 0s | < a | Neither can it 
be eliminated by because S* cannot be sustained for a long time and f is a long- 
term force. From another point of view, since 0 a is simulated by 0s and g^ and g® are 
just g^ rotated by and P^ respectively, it is reasonable to simulate g* by g^ . 

The block diagram for the revised algorithm is shown in Figure 5.2. g'^ is no 

longer followed but g®^ is used as direct input, g®*^ is g in an imaginary reference 
frame FrsR. If the simulator only responds to angular input, the simulator frame is FrsR. 

SR S 

In the pure rotational case, FrsR= Frs, then g = g 

The original algorithms had an active translational and tilt channel under pure 

rotational input. Neither P^^ nor was zero and it was often unstable imder pure 

rotational input or a mixture of translational and rotational input The revised algorithm 
generates the same results as the results generated by the original adaptive algorithms 
under pure translational input. The difference is that the revised algorithm has a null 

translational channel under pure rotational input. Both ^ and ^ are zero. The revised 
algorithm is stable under translational, rotational, or a mixture of translational and 
rotational inputs. 
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Desired to be in Frs The simulator will try to compensate 

i I I Fr$ it by S*’ and simulator tilt 



Figure 5.2. Revised UTIAS Adaptive Algorithm. 














5.2. Spikes in the Outputs of the NASA and TJTIAS Adaptive Algorithms 

Some significant spikes occurred when either the NASA Adaptive algorithm or 
the UTIAS adaptive algorithm was run as shown in Figure 5.4. Whenever there was a 
sharp change in the translational acceleration input, a spike would occur. In the UTIAS 
report [1], [2], the spikes were said to be due to the attenuation of the second pulse of a* 
(simulator acceleration). But after careful examination, the real cause of the spikes is 
identified to be the rotational acceleration of the simulator. Since the pilot’s head is in a 
position some distance away from the centroid of the simulator (the rotational center in 
the adaptive algorithms), a simulator rotational acceleration will generate some additional 
acceleration at the pilot's head. It is this acceleration that caused the spikes. 

The specific force on the pilot’s head can be expressed as; 

( 5 - 2 - 1 ) 

where a^ is the acceleration at the pilot's head, a^ is then expressed as 

§re= a" + Es + 2 ffls X + Q® xRg + ©j x(fflj X Rj,) (5-2-2) 

where Rs is the vector from the centroid of the simulator to the pilot's head. 

Assuming = Rg = 0 , then 

3re= + ®s X Es + ®s ^ (®s X Es) (5-2-3) 

Equation (5-2-3) indicates that both © ® and © g contribute to a|g . The limits for tilting 
^ are set to 0.0524 rad/sec. Usually each component of ^ is less than 2 m for most 
simulators. Therefore each component of ©g x (oo® x Rg) has a magnitude less than 
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0.05 m/sec^. This tOTn will not contribute much to . Major attention will be paid 
to fflg , which is expressed as; 


(pi + qj + rk)x(Rs^i + Rsyj + Rs^k) 

= (Rs,q-Rsyr)i + (Rs*r-RszP)j + (RsyP-Rsx^k 


(5-2-4) 


There are two ways to prove that x ^ is the cause for those spikes. First, the 

magnitude of o g ^ Es estimated by hand calculation. Then the value can be 

compared with the magnitude of the spike. The second approach is to set R = (0,0,0), 
when the spikes are supposed to be eliminated completely. These two tests were 

performed and proved that ® g x Rg is the real cause of the spikes. 

Since the spikes are generated by ^ x ^ , there are two ways to decrease the 

spikes. The first is to reduce ® g and the second is to reduce Rg. ^ is determined by the 
position of the pilot's head relative to the position of the center of the simulator rotation. 
Only the former can be reduced. In the UTIAS report [3], [4] the center of rotation is 
selected to overlap with the centroid of the simulator for minimum extensions of actuator 
1^. This selection is quite reasonable and any change will increase the actuator 

extensions. Then the only approach left is to reduce 6 g . 

The generation of was investigated as shown in Figure 5.2. The Pitch/Surge 

filter was cammed as an example. The Roll/Sway filter has similar characteristics. In 
the Pitch/Surge filter, 

q =e, =px2*4+px3*e. (5-2-5) 
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A rotational acceleration will be generated by: 

q - e. = • (fL(i+l) - £L,(0 ) / * . (5-2-6) 

where i indicates the current time step and i+1 the next time step. If there is a sudden 
change in an extremely large 0, might be generated. Furthermore, the smaller the 
step size dt is chosen, the larger 0, will be generated. For example, if the longitudinal 
acceleration input is a step or a pulse input, and when t = 0, 0, =0, then at the first time 
step 0j = px2 * dt. The UTIAS adaptive algorithm set p *2 = 0.12 initially and dt = 
0.05 sec. Then an input which has an amplitude of 1 m/sec^ would generate 0^ = 0.12 * 
1 / 0.05 = 2.4 rad/sec^ at the first time step. The UTIAS simulator used R = -0.2 i - 0.465 
j - 1.783 k. Then q = 0, = 2.4rad/sec^ would generate: 

® s ^ Es = (Rszq - RsyO 1 + (Rsxf - RszP) J + (RsyP " ^ 

= Rs^qi-Rs,,qk 

= -1.783-2.4 1 + 0.2 -2.4 k 

w - 4.3 i + 0.48 k (m/sec^) 

© g X Rg contains a component with a magnitude of about 4 m/sec^ in the x direction at 
the pilot's head. This acceleration will overpower all other simulator motion effects and 
give the pilot a perception in the wrong direction at the begirming of the simulation. This 
significantly wrong perception will also happen whenever there is a sudden change in 
aircraft translational acceleration. 

A limit on tilting angular acceleration © of 0.5 rad/sec^ has been added in the 
current software. Spikes were attenuated significantly at the price of a slower simulator 
tilt response as shown in Figure 5.5. In the former example, if the limit of 0, is set to be 
0.2 rad/sec^, then the magnitude of the spikes can then be attenuated by about ten times. 
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At the same time, about 10 * dt is needed to get 0^ to reach the expected value. Then at 
about 0.5 seconds, the simulator response is slower than before the new limit is added. 
Since it only influences about 0.5 seconds duration, this is more tolerable than the spikes, 
which completely reversed the pilot's perception. On the other hand, simulator 
translational motion performs a more important role at the banning of simulating a 
translational acceleration change. This makes the slower response more tolerable. 
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Atrcraft Specific Force at Pilot Heaci 




Aircraft Seneed Spedie Force at PiM Heed 



Figure 5.4. Specific Forces Generated by the NASA Adaptive Algorithm with No Limit 
on Angular Acceleration. 
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Aifcmlt Spedfic Force et Pilot Head 




Figure 5.5. Specific Forces Generated by the NASA Adaptive Algorithm With 0.5 
rad/sec/sec Limit on Angular Acceleration. 
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6. Phase 1 Results 


6.1. Translational Input Case 

Both a pulse input of 1 m/sec^ for a 1 0 second duration and a sinusoidal input of 1 
m/sec^ with a frequency of 3 rad/sec were used as translational inputs. In each case the 
input was applied individually to each translational degree of freedom (surge, sway, and 
heave). The pulse contains both high frequency and low frequency components. Figure 
6.1 shows the specific force at the pilot’s head for the pulse input, and Figure 6.2 shows 
the response for the sinusoidal input. 

The response onsets, i.e. the start of a response, generated by each algorithm are 
significantly different During the response-sustained period the difference among those 
results is not significant. Major attention will be given to the response onsets. For 
comparison of the response onsets Figures 6.1 through 6.3 were generated with a shorter 
time axis (5 seconds). 

As observed in Figure 6.1 the classical algorithm and the UTIAS adaptive 
algorithm response onsets have smaller magnitudes and larger phase lag as compared to 
the NASA adaptive and optimal algorithm response onsets. If the responses within the 
first second are inspected, the classical algorithm generates a response with a magnitude 
equal to about one half the magnitude of the aircrafi response. The other algorithms 
generate responses with higher magnitude. Note that for an easy comparison the input 
scaling gain has been set to one for all algorithms. The optimal algorithm had the least 
onset lag and the NASA ad^tive algorithm had the second least onset lag, with the 
classical algorithm and the UTIAS adaptive algorithm having larger lags by about 0.1 to 
0.2 seconds. The “sags” following the onsets generated by the classical algorithm and the 
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UTIAS adaptive algorithm were significantly larger than the sags generated by the 
optimal and the NASA adaptive algorithm. Large sags usually decre^e the fidelity of 
motion cues. 

The response offsets are similar to the response onsets. It should be noted that the 
UTIAS adaptive algorithm has some excessive delay. This delay is due to the tilt 
ad^tive parameter p 2 being driven to a very small value after the input magnitude is 
sustained for several seconds. The NASA adaptive algorithm has a fixed tilt channel 
parameter ks so that it avoids the extra delay. 

These differences on response onsets and offsets might have significant meaning 
in the overall fidelity of those algorithms. For further comparison, in Figure 6.2 the 
specific force at the simulator pilot's head corresponrUng to sinusoidal input was 
generated. It is observed that the NASA adaptive algorithm and the optimal algorithm 
generated responses with high magnitude. The UTIAS adaptive algorithm generated 
response with medium magnitude and the classical algorithm generated response with 
low magnitude. The optimal algorithm has the highest fidelity in shape. All other 
algorithms generated some extra extremum points. The NASA adaptive and UTIAS 
adaptive algorithms generated specific forces in a wrong direction during about the 
beginning 0.3 seconds. This problem on specific force in a wrong direction was 
discussed in Section 5.1. 

On the other hand, the optimal algorithm and NASA adaptive algorithm resulted 
in significantly larger actuator extension costs than the classical algorithm and the UTIAS 
adaptive algorithm. 
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A/C & Simu. Specific Force atPilotHead 




Line with no marks -- Aircraft response. 

Line with marks - Simulator response. 

Classical o~ NASA Adaptive +-- UTIAS Adaptive x- Optimal 

Input: pulse x-acceleration input, magnitude =1.0 m/sec^, duration =10 seconds. 

Figure 6.1. Specific Force at Pilot's Head Corresponding to Pulse Input. 
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A/C & Simu. Specific Force at Pilot Head 
1.5 I , — . 



t (sec) 


Sensed Specific Force 



Line with no marks -- Aircraft response. 

Line with marks — Simulator response. 

Classical o~ NASA Adaptive +-- UTIAS Adaptive x— Optimal 

Input: sinusoidal x-acceleration input, magnitude =1.0 m/sec^ © = 3 rad/sec. 

Figure 6,2. Specific Force at Pilot's Head Corresponding to Sinusoidal Input. 
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6.2. Rotational Input Case 

A doublet pulse angular acceleration input of 0.05 rad/sec^ for a 10 second 
duration was used. The classical angular rate response has some overshoot when it drops. 
This causes the simulator angular position response to have a slight imwanted drop 
(Appendix B Figure B.l). The amplitude of the angular position generated by the 
optimal algorithm slightly exceeds the aircraft angular position (Appendix B Figure B.2). 
The UnAS adaptive algorithm and the NASA adaptive algorithm both have smooth 
angular position response. The magnitudes of the simulator angular position responses 
are significantly attenuated, especially in the NASA adaptive algorithm (Appendix B 
Figures B.3 and B.4). In the UTIAS Report Volume 2 [4], the Bode fi^equency responses 
of the classical algorithm and optimal algorithm were presented. These plots show that 
the gain of the classical and optimal respoiLse fluctuates around one and some small phase 
shift would happen in both the classical algorithm and optimal algorithm. These analyses 
are consistent with the graphs in Appendix B. For ftirther comparison some sinusoidal 
input cases were studied. From Figure 6.3 it can be seen that in some input cases the gain 
of the classical and optimal algorithms could exceed one. The phase shifts of the 
classical algorithm and optimal algorithm are larger than the phase shifts of the UTTAS 
adaptive algorithm and NASA adaptive algorithm. The NASA adaptive algorithm has a 
response with much smaller magnitude. 

The optimal algorithm resulted in file largest jack extension costs. The UTIAS 
adaptive algorithm and the NASA adaptive algorithm resulted in jack extension costs 
smaller than the classical algorithm and the optimal algorithm. 
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A/C and Simu. Sensed Angular Rate 



t (sec) 

Line with no marks — Aircraft response. 

Line with marks — Simulator response. 

Classical o— NASA Adaptive +— UTI AS Adaptive x— Optimal 

Input: Sinusoidal roll acceleration input, magnitude = 0.05 rad/sec^, <0 = 3 rad/sec. 

Figure 63. Aircraft and Simulator Angular Rate Response. 
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6.3. Conclusions and Recommendations 


At the current stage, in translational input cases the optimal algorithm and the 
NASA adaptive algorithm perform better than the classical algorithm and the UTIAS 
adaptive algorithm. In rotational input cases, the UTIAS adaptive and the NASA 
adaptive algorithm perform better than the classical and the optimal algorithm. 

The classical algorithm is a type of linear algorithm distinguished by the fact that 
the washout output is proportional to the input. This is an imdesitable attribute because 
the duration of the onset cue is then limited by the maximum amplitude input case. 
Therefore in lower amplitude case the onset cue cannot be sustained any longer even 
though the hardware motion resource is quite available. The classical algorithm is 
therefore not recommended. 

The NASA adaptive algorithm has been improved by eliminating imwanted 
spikes and excessive oscillations. The NASA adaptive algorithm currently performs 
well. 

The UTIAS adaptive algorithm is expected to perform better than it does 
currently. It has also been improved by eliminating imwanted spikes and excessive 
oscillations. Some more improvements are needed such as increasing the slow tilt 
response by fixing or restricting the tilt adaptive parameters, thus changing some 
parameters to make the translational response be sustained for a longer time. With these 
further improvements implemented, the UTIAS adaptive algorithm might be attractive. 

For the optimal algorithm, if selecting another washout frame can decrease the 
actuator extensions, this algorithm will be very attractive. Currently the optimal 


75 



algorithm performs washout in the pilot head frame. As discussed in the UTIAS report, 
this frame selection results in some extra jack extension costs when the simulator rotates. 

A washout frame whose origin lies on the simulator centroid is desired. But the 
desired frame will cause some cross coupling in the pilot specific force sensation. This 
cross coupling might make the generation of the optimal transfer frmctions much more 
complicated or even not practically feasible. Further study on the optimal washout in the 
simulator centroid frame is recommended. 
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7. Phase 2 Optimal Algorithm Redevelopment 

In Phase 2, further comparisons are made between the NASA adaptive algorithm 
and a redeveloped optimal algorithm. The classical algorithm will not be evaluated due 
to the undesirable results obtained in Phase 1. Since the adaptive algorithms are 
fundamentally similar with the NASA adaptive algorithm having more desirable features 
than the UTIAS algorithm, the UTIAS adaptive algorithm will not be evaluated in Phase 
2 . 

The optimal algorithm will be revised with the center of rotation redefined and 
utilizing a new cost function with additional terms. In addition, a new design approach 
that integrates FORTRAN/MATLAB/SIMULINK is developed in which a set of selected 
parameters can be tested in 30 seconds as compared to 15 minutes with a conventional 
approach. The revised optimal algorithm also includes a new vestibular model. 

Two general problems that were unresolved in Phase 1 are addressed in Phase 2. 
The first problem was that for inputs with large magnitudes all the algorithms tended to 
drive the simulator beyond its motion limit. A nonlinear gain algorithm was developed 
so that the simulator would not reach its motion limit for most input cases. The second 
problem was that when the simulator reached its limit an algorithm was needed to brake 
the simulator to a full stop and release the brake at a proper time to allow the simulator to 
again follow the washout algorithm output. For this problem a braking algorithm was 
developed in Phase 2. 
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7.1. New Vestibular Models 


The semicircular canals have been well studied by many researchers. A model 
that was more consistent with the results of recent research could be found through a 
review of the reports written by those rwearchers. The construction and justification of 
this proposed model is discussed in detail in Appendix A. The revised model for 
practical usage (neglecting the short time constant) can be expressed as: 

o T^T,5(1+T,:y) 

' ® (1+T,S)(1+T^5) ^ 


where o is the input angular acceleration and ® is the sensed angular velocity. x„ tl, 
and Xp are time constants. Xa = 30 sec, Xl = 0.06 sec, and Xp = 5.7 sec. 

The UTIAS report [3] suggested a simplified otolith sensation model that ignored 
the short time constant term (XgS + 1) in Equation (2-1-2) for practical usage: 


f _ K(x3S + 1) 

i (TlS + I) 


(7-1-2) 


A 

where f is the input specific force and f is the sensed specific force. 

Human sensation models (7-1-1) and (7-1-2) are used in the revised development 
of the optimal algorithm. 


78 



7.2. Redefinition of Center of Rotation 


In Phase 1, the optimal algorithm resulted in significantly, larger actuator 
extensions than the other algorithms. It was found that this problem arose from the 
different selection of the center of rotation. If the center of rotation is selected at the 
centroid of the simulator, there are several cross-couplings; 

pitch surge, pitch heave, roll sway, roll heave, yaw -> sway, yaw surge. 
This can be shown in Figure 7. 1 : 


y(a?,R,) 



Ct (Centroid of the simulator) 

Figure 7.1. Phase 1 Optimal Algorithm Center of Rotation 
An angular velocity ca g with respect to the centroid of the simulator always generates a 
translational motion at the pilot’s head. 

If the pilot's head is selected as the center of rotation and the vestibular system is 
considered as the only motion sensation organ of the pilot, all the cross-couplings will 
disappear. This is the reason for which the UTIAS report [3] chose the pilot’s head as the 
center of simulator rotation. Unfortunately, this selection resulted in excessively large 
actuator extensions in some input cases. This can be shown by an example in which the 
simulator has both translational and rotational displacement.. For example, in the 
pitch/surge mode the same tilt angle 0 and the same linear displacement S* are generated 
in two different systems, one of which has the pilot’s head and the other has the centroid 
of the simulator as the center of rotation. The motion of the simulator is shown in Figure 
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7.2. Figure 7.2 shows clearly that the different selections of center of rotation will result 
in different amounts of actuator extensions. 

The adaptive algorithms selected the centroid of the simulator as the center of 
rotation. However, cross couplings were ignored in the algorithm development. Test 
results showed that some undesirable spikes were generated due to the cross couplings. 
Some effort was spent on reducing or eliminating the magnitudes of the spikes in Phase 
1. The spik« were reduced but not eliminated. The spikes may generate some negative 
motion cue to the pilot’s sensation or at least cause some delay on the pilot’s sensation on 
the positive motion cue. In Phase 2, a new optimal algorithm that chooses the centroid of 
the simulator as the center of rotation was developed. Cross couplings were handled 
explicitly and the spikes were eliminated in the new optimal algorithm. 
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Centroid of simiilator is the center of rotation 


Pilot’s head is the center of rotation 
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the original position 

E] and E3 have the centroid of the 
simulator as the center of rotation 

E2 and E4 have the pilot's head as the 
center of rotation 
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7.3. Revised Development of Optimal Algorithm 

7.3.1. Construction of the System Equations and Cost Function 

The optimal filters for the four modes, pitch/surge, roll/sway, heave, and yaw are 
designed separately in the optimal algorithm development. This design technique was 
also employed in the development of the other three motion-base drive algorithms. 

In Phase 2 only the pitch/surge mode is redefined. The roll/sway mode filter 
design was based on the pitch/surge redevelopment. The heave and yaw modes are 
unchanged from Phase 1. The revised development for the pitch/surge mode is as 
follows: 

(a) Rotational Motion Sensation 
Input: 


u = 


1 

< D \ 
1 


1 

1 

k. 


“ 2 . 


The semicircular canal sensation model given in Equation (7-1-1) is then used to obtain 
the rotational sensation q : 


q = 


T,S (1 + Ti^s) 


(1 + T3SXI + XpS) + TjS + To 


TjS^ + TjS 


U, 


where T„ = T, = = TlX^Tq, T, = x.Tq 


(7-3-2) 


■Cl^a 




Now expressing the above in state space notation yields the following equations: 
Define Xi -2 = [^1 ^ 2 ]^’ fhen 


|x,-2 = A„,x,.2+B„3U 

1 q=C.„x,_2+D„,u 


(7-3-3) 
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where 


-T, 

1' 

, B = 

X-TJ, 

O' 

-To 

0 

9 ^SSC 

-t„t. 

0 


= 

C„=[l o]. d,.=[t, «] 

(b) Translational Motion Sensation 

The otolith sensation model given in Equation (2-1-2) can be rewritten in terms of the 
break frequencies Ao and Bo; 


f 


Gq(s + Aq) 

(s + Bq) 


(7-3-4) 


where f is the input specific force and f is the sensed specific force. Go - 5.86, 
Ao = 13.2 sec. Bo = 5.33 sec . 

Input ; 


u = 


1 

<D: 
1 


1 

LaJ 


1 

c 

N) 

1 


The center of simulator rotation is set at the centroid of the simulator. Then 

f = a,+g0-Rs,0 (7-3-5) 

where Rsz is the z-component of the vector from the centroid of the simulator to the 
pilot’s head. Equation (7-3-5) can also be expressed in the LaPlace domain; 


Specific force sensation: 


f = Go^^f 
°s + Bo 

S “h A 1 


- G.[- 


- RszS^ - RszAqS^ + gs + gAp s + At 


s^+Bps" 


S + Br 


]U 


^ p r Rsz(Bo-Ao)s"+gS+gAo 

s' + Bps" 


State-space model: 

Define = [x, x^ Xj x^]^, then 


Ap - Bp 


(7-3-7) 


= A.^x^6 + B^u 

I f “ ^oto?3~6 '*' ®otoil 


(7-3-8) 


where 


A_ = 


‘oto 


-Bp 

0 

1 

0 


Bj* (Bq Ap) 

0 

0 

0 

0 

0 

, B^ = G„ 

gAp 

0 

0 

1 

0 

0 

’ oto 0 

g 

0 

0 

0 

0 

-Bo. 


0 

Ap — Bp 


C„^=[l 0 0 11, D^=Gp[-r. 1] 

(c) Combination of Both Rotational and Linear Motion Sensation 


JXk6 = AvX,~6 + BvU 
[ y, =CvX,-6 + DvU 


(7-3-9) 


84 



where 


Xl ~6 = [^1 ^2 ^3 X4 X5 Xj , Yj = 


q 

yi. 


Av = 


5SC 

0 A. 


B. 


B. 

B. 


Cv = 


ISC 

0 c 


0 


, Dv = 


oto_ 


D. 

D 


oto_ 


It is assumed that the same human motion sensation model can be applied to both 
the pilot in the aircraft and the pilot in the simulator. Define the input to the pilot in the 

aircraft as u^ = 


0 . 


u 


A 1 


Ax . 


U 


and the input to the pilot in the simulator as 


A 2 . 


U, = 


e. 


■'Sx 


U 


SI 


U 


S 2 


To reduce the order of the system equation, defme x« = “ 2 ^a > 


and then define the pilot’s sensation error as e = y, - v , then 

— IS — lA 


*e ~ ~ AyXg+BvUg (AyXA+ByU^) 

= Ay X^ + By Ug — ByU^ 


( 7 - 3 - 10 ) 


e=CyXe +ByUg -DyU^ ( 7 - 3 - 11 ) 

(d) It is necessary for the control algorithm to explicitly access some motion states, such 
as the linear velocity and displacement of the simulator, which are desired to appear in 
the cost function. For this purpose some additional terms were grouped and included in 
the system equations. 




'x/ 



Xg 

Ja,dt 

— 

X9 

0 


X,0 

0 


-Xii_ 


( 7 - 3 - 12 ) 
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2£d ®d ii 


where 


"0 

1 

0 

0 

o' 


"0 

o' 

0 

0 

1 

0 

0 


0 

0 

0 

0 

0 

0 

0 

> B,= 

0 

1 

0 

0 

0 

0 

1 


0 

0 

0 

0 

0 

0 

0_ 


_1 

0 


(e) Input u consists of filtered white noise. This can be expressed by the state-space 
equation 


= A„ x„ +B„ w 

1 M = x„ 

where 


A. = 

■-1.0 

0 ■ 

, B = 

1.0' 

D 

0 

-1.0 

’ O 

1.0 


( 7 - 3 - 13 ) 


(f) Assemble system equations. 

The desired system equation can be obtained; 
Define 


X = 
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then 


X = Ax + Bu, +Hn 
e 


y = 


5dj 


= Cx + Du, 


where 


1 

> 

< 

0 

1 

< 


1 

W 
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0 A, 0 

, B = 


, H = 

0 

1 

O 

o 

> 

9 

L 


1 

o 

1 


r 

a 

s 

1 


C = 




with the cost function 


(7-3-14) 


J = E|{‘(e^Qe-HxlRaXd +ujRus)dt} (7-3-15) 

The cost function implies that three variables are to be constrmned; e, Xd, and us. 
e is the sensation error, xa, which is a group of terms, along with us define the rotational 
and linear motion of the simulator. The cost function constrains both the pilot sensation 
error and the simulator motion. 

7.3.2. Transformation to Standard Form 

The system equations and cost function can be transformed to the standard form 
by the following equations: 

System equation: x - A'x + Biri -i- Hn (7-3-16) 

Cost function: J= E|£[x’^R;x + u:%u:]dt| (/J=0) (7-3-17) 
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where 


G = 


, Ri =C^GC, R„=C^GD, Rj =R + D^GD, 


Q 0 

0 R, 

A' = A - BR-'Rjj , u: = u, + Rj'R^iX , R| = R, - RjjRj’R^ , 


u: = -R-' P(t) X 


7.3.3. Solution of the Optimal Control Problem 

The standard form solution to the optimal control problems can be applied since 
both the system equations and the cost function are in standard form. Therefore, 

(7-3-18) 

where P is the solution of the algebraic Riccati equation; 

0 = R; -PBRj' B^P + A'^'P + PA' (7-3-19) 

From (7-3-17) and the definition of u^, it is obvious that 

Us = -[r;‘(B^P + R,^,)]x (7-3-20) 

Define 

F = R-'(B^P + Rjj) (7-3-21) 

Then 


Us = -Fx 


Partition F corresponding to the partition of x = 







L^n. 



1 

IX 

rt 

1 


J 

X 

u.=-[f. F, F,] 

Xd 

II 

1 




1 

c 

XI 



d _ 


(7-3-22) 


(7-3-23) 
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Eliminate from the system equations; 









'Ay 0 

-By- 

c 

Xh 

+ 

By 

_ 0 A, 

0 

— a 




(7-3-24) 


After taking the Laplace transform of (7-3-23) and (7-3-24) and making some 
substitutions, the following equations are obtained: 


u,(s) = W(s)u.(s) 


(7-3-25) 


where 


W(s) = [f, 


^Sl- Ay +ByF, 

ByF, 

-1 

By(I + F 3 ) 

{ b,f, 

sI-A^ +BjFj 


B-B, 


Note that the filter matrix W(s) is not identical to the filter matrix W(s) described 

in Section 4 and implemented in the UTIAS report. The revised design of W(s) was 
based on control inputs of angular acceleration and linear acceleration, while the design 
of W(s) was based on control inputs of angular displacement and linear acceleration. A 
revised on-line Fortran development was not implemented; the pitch/surge and roll/sway 

filters W(s) developed in Section 4 were simply replaced with the revised filters W(s) in 
the Phase 1 on-line development. 

7.3.4. Solution of the Washout Filter Coefficient 

W(s) is solved by Matlab functions and Fortran programs. All the system 
matrices were constructed in Matlab functions. First, the weighing matrices Q, R, and Rd 
in the cost function were selected. These matrices were adjusted several times in a trial 
and error procedure. Then the system matrices Ay, By, Cy, Dy, Ri, R2> Ri2» A', and R[ 
were constructed in a proper order. The algebraic Riccati equation was defined in the 
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Matlab ftinction. The definition was transported to a Fortran environment in which the 
algebraic Riccati equation would be solved by Fortran programs. The solution of the 
algebraic Riccati equation was transported back to the Matlab environment and Fi, Fi, 
and F 3 could then be defined. 

After all matrices in (7-3-25) were defined, the solution of W(s) still required 
symbolic calculation since there was the Laplace variable s. A series of fairly 

complicated Fortran programs were presented in the UTIAS report to solve W(s). The 
solution required manual cancellation of common poles and zeros. It was found that the 
problem could be solved more efficiently by Matlab after a careful observation of 
Equation (7-3-25). Recall that a linear control system can be expressed in state-space 
form: 


Jx = Ax + Bu 
|y = Cx + Du 


(7-3-27) 


where A, B, C, D are system matrices. 

Equation (7-3-27) determines a transfer function between y and u: 

y = [C(sl - A)"‘ B + D] u (7-3-28) 

By an observation of the similarity between (7-3-26) and (7-3-28), it was found that the 
state-space model that had ua as the input and us as the ou^Jut could be obtained directly; 


[ X = Ax + Bu^ 
Us =Cx + Du^ 


(7-3-29) 
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\^iiere 



-BvF 

a.-b,f 


, B = 


Bv(I + F3> 


C = [F, Fj, D = -F3 


This state-space model could be converted to the transfer function model by the Matlab 
function ss2tf easily. A Matlab function was also written to cancel common poles and 
zeros and yield some transfer functions in reduced order. A simulation setup that reads in 

the transfer functions automatically was constructed in Simulink. The effects of W(s) 

are visualized with the simulation. If W(s) is not satisfactory, the procedure of the 

design of W(s) is repeated by selecting some new cost function matrices Q, R, and Rd, 
until a satisfactory result can be finally approached. The procedure for the solution of 

W(s) is illustrated in Figure 7.3. 
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Ouq>ut optimal filters 
to tfiis.dat 


Figure 7.3. W(s) Solution Procedure. 
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7.4. Nonlinear Gain Algorithm 

It is desirable to maximize the magnitude of the motion cue while remaining 
within the operational limits of the motion system. Although it is very difficult to restrain 
the response within hardware limits all the time, it was found feasible to restrain the 
response successfully during most of the simulation time by implementing a nonlinear 
gain in the simulator motion-base driving software. 

When the magnitude of input to the simulator motion system is small, the gain is 
desired to be relatively high, otherwise the output may be below the pilot’s perception 
threshold when it should be perceptible. When the magnitude of input is high, the gain is 
desired to be relatively low otherwise the simulator may go beyond the hardware limits. 
Define the input as x and the output as y. Define Xmax as the expected maximum input 
and ynux as the maximum output, and So and Si as the slopes at x = 0 and x = Xnux 
respectively. Four desired characteristics for the nonlinear gain can be expressed as; 

(1) X = 0 => y = 0 ; 

(2) x=x^ => y= Ym® ; 

(3) y’L”S„; 

w y = s, ; 

A third order polynomial can be employed to provide functions with all the desired 
characteristics. This polynomial will be of the form; 

y = CjX^ + CjX^ + c,x + Co 


93 



where 


Co = 0 

C, =So 

C 2 = X™ • (3y™„ - 2SoX^ - S,x^ ) 

C3 = • (SqX^ - 2 y^ + S,x^ ) 


One example of this nonlinear gain is shown in Figure 7.4, with parameters set as Xn»x = 

10, ym« = 6, so= 1.0, Si = 0.1. 


Nonlinear Gab) 



Figure 7.4. Nonlinear Gain. 


The nonlinear gain developed by NASA for its motion-base drive software results 
in a high probability of reaching system limits for accelerations higher than 3 m/sec^, 
which is common in aircraft simulation. When the new nonlinear gain as shown in 
Figure 7.5 was implemented in the software, in most input cases the simulator motion 
was within the hardware limit. 
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7.5. Actuator Extension Limiting 

There is a possibility that the output of the washout algorithm will drive the 
simulator beyond its hardware limit. It is obvious that any hardware shut down due to the 
simulator’s excessive excursion should be avoid^. This raises the requirement that the 
simulation software should be able to handle the situation that the simulator may need to 
be arrested before encountering the system limit. This necessary software subroutine is 
called the braking algorithm. When necessary, the braking algorithm takes over the 
control of the simulator from the washout algorithm. It is also desirable that when the 
washout algorithm would begin to drive the simulator toward smaller excursions, the 
braking algorithm returns control to the washout algorithm. In other words, the braking 
algorithm should release the system to resume r^lar simulation at the proper time. The 
logic of the braking algorithm is shown in Figure 7.6. 

The algorithm makes a series of decisions to determine when to brake the 
simulator and while braked to determine when to release the brake. The first decision is 

based on an evaluation of the expression (2-Co-ab-s — v^) for each actuator at each 
simulation cycle, where v is the velocity of the actuator, s is the available actuator stroke, 
3b is the acceleration/deceleration by which the braking algorithm will stop the actuator, 
and Co is a coefficient less than or equal to 1 that is described below. 

3 ), should not be larger than the maximum actuator acceleration ^ich is a 
hardware parameter. The simulation software sets a software limit position that is before 
the simulator hardware limit position. If the simulation cycling rate is infinitely high, 
when Co = 1, the actuator will be stopped exactly at the software limit position; when 
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Co < 1, the actuator will be stopped before the software limit position. Since the real 
system always has a finite cycling rate, it is possible that the actuator may go a little 
beyond the software limit position even though a braking algorithm is implemented. For 
this reason, the software limit position should be set slightly before the hardware limit 

position. When one actuator reaches the braking region defined by(2-Coab-s-v^<0) 
at time to, it will be decelerated by ab. At the same time, all other actuators will be 
decelerated proportionally to their respective velocities at time to. 

The second decision is based on a comparison between the actual simulator sta te s 
and the washout algorithm output. The simulator states refer to the simulator’s linear and 
angular positions. When the simulator has been completely stopped by the braking 
algorithm, the comparison begins. When the washout algorithm output states are smaller 
than the corresponding actual simulator states, the simulator will begin to follow the 
washout algorithm output again. 

When the braking algorithm releases the brake, the washout algorithm output may 
have large velocities while the simulator has small velocities at that instance. To avoid 
the regular simulation being resumed with any excessively large simulator acceleration, 
an algorithm which allows the simulator to follow the washout algorithm output 
gradually was designed. The algorithm expresses the linear and angular position 
command to the simulator as: 
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Xc_ = »(t) • Xa. + (1 - m ■ (Xw.i - Xa„) 


( 7 - 5 - 1 ) 


where 


2£com the commanded position which drives the simulator; 
X^ct current actual position of the simulator; 

2^w«di current output of the washout algorithm. 

a(t) is a function of time. It increases from 0 toTi / 2 gradually. 
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n>=0 INDICAIXS THAT THE OUTPUT OF THE 
WASHOUT FETER IS TOWARDS SMALLER BETAS 
AND SSI. ID IS DETERMINED IN THE MAIN ROUTTNE. 


YES 

j RETURN RL(I) 


Figure 7.6. Braking Algorithm in the Actuator Driving Subroutine. 
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8. Phase 2 Results 


8.1. t>itch/Surge Mode and Roll/Sway Mode 

The pitch/surge mode and the roll/sway mode have similar characteristics with 
both employing the simulator tilt to provide sustained specific force cue, and both having 
the cross-coupling from the simulator tilt to the pilot’s head translational motion. Both 
half cycle sinusoidal inputs and ramp to step inputs were employed to test the pitch/surge 
and roll/sway channel of the NASA adaptive washout algorithm and the optimal washout 
algorithm. 

When a 0.1 Hz half cycle sinusoidal translational acceleration was employed as 
the input, the output of the NASA adaptive algorithm and the output of the optimal 
algorithm were similar. The specific forces at the simulator pilot’s head generated by the 
two algorithms had similar magnitudes and shapes. The actuator extensions resulted 
from the two algorithms were also similar. One difference with the NASA adaptive 
algorithm was that the specific force generated went to a wrong direction at the beginning 
of the input and had an extra hump besides the sinusoidal hump corresponding to the half 
cycle sinusoidal hump of the input. The specific force output of the optimal algorithm 
did not have these two deviations from the shape of the input. The outputs are shown in 
Appendix C, Figures C.l, C.2, and C.3. 

When a 0.5 Hz half cycle sinusoidal translational acceleration was used as input, 
the shape of the specific force at the simulator pilot’s head generated by the NASA 
adaptive algorithm had some significant distortion. The specific force in the wrong 
direction at the beginning of the input had a magnitude which was about half of the 
magnitude of the onset of the specific force cue and a duration which was about one half 
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of the time by which the input reached its maximum. The specific force onset had some 
more wiggles than the input and the washout of the specific force had a large ova^hoot. 
Stimulated by the same input, the optimal algorithm generated a specific force at the 
simulator pilot’s head with significantly smaller shape distortion. There was no negative 
cue at the beginning of the input The overall shape of the output was similar to the 
input, with a small washout overshoot The outputs are shown in Appendix C, Figures 
C.IO, C.ll, andC.12. 

When a ramp to step acceleration with a ramp slope of 3 m/sec^ was employed as 
the input, the specific forces at the simulator pilot’s head generated by the two algorithms 
had similar magnitudes and shapes. The actuator extensions resulting from the two 
algorithms were also similar. The specific force generated by the NASA adaptive 
algorithm went in the wrong direction at the beginning of the input The specific force 
generated by the optimal algorithm reached its ma ximum earlier as compared to the 
adaptive algorithm and then slowly decreased. The outputs are shown in Appendix C, 
Figures C.4, C.5, and C.6. 

When a ramp to step acceleration with a ramp slope of 5 m/sec^ was employed as 
the input, the shape of the specific force at the pilot’s head generated by the NASA 
adq>tive algorithm had some significant distortion. The specific force in the wrong 
direction at the beginning of the input had a magnitude about one half of the magnitude 
of the onset of the specific force cue with its duration almost equal to the input ramp 
duration. The onset of the specific force was followed by a sag. Stimulated by the same 
input, the optimal algorithm generated specific force at the pilot’s head with significantly 
smaller shape distortion. There was no negative cue at the beginning of the input. The 
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overall shape of the output wm similar to the input. The ou^uts are shown in Appendix 
C, Figures C.7, C.8, and C.9. 

8.2. Heave Mode 

A pulse acceleration with a magnitude of 3 m/sec^ and duration of 5 seconds was 
employed as the input. The onset of the specific force at the simulator pilot’s head 
generated by the NASA adaptive algorithm was similar to the one generated by the 
optimal algorithm. When the specific force at the airo-aft pilot’s head decreased, the 
optimal algorithm generated a corresponding drop of the specific force at the simulator 
pilot’s head, while the NASA adaptive algorithm did not generate an obvious drop of the 
corresponding output. The drop of the specific force at the simulator pilot’s head 
generated by the optimal algorithm was followed by a washout with a large overshoot. 
The outputs are shown in Appendix C, Figures C.13, C.14, and C.15. 

8.3. Pitch Mode and Roll Mode 

A doublet angular acceleration with a magnitude of 0. 1 rad/sec^ and a duration of 
5 seconds was employed as the input to the pitch channel and roll chaimel. For the pitch 
test runs the shape of the specific force generated by the NASA adaptive algorithm was 
nearly the same as the shape of the input, while the specific force generated by the 
optimal algorithm slowly decreased. For the roll test nm the specific force generated by 
the NASA adaptive algorithm was nearly the same as the shq)e of the input, while the 
specific force generated by the optimal algorithm required more time to settle during 
washout. Some distortion with the simulator angular velocity relative to the input was 
also observed with the optimal algorithm. The outputs are shown in Appendix C, Figures 
C.16, C.17, C.18, C.19, C.20, and C.21. 
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8.4. Yaw Mode 


A pulse angular acceleration with a magnitude of 0.05 rad/sec^ and duration of 5 
seconds was employed as the input. Both the outputs generated by the NASA adaptive 
algorithm and the optimal algorithm have the characteristic that the simulator yaw rate 
and yaw angle were washed out toward zero when the aircraft yaw rate was rising up to a 
constant value and the aircraft yaw angle was always increasing. The simulator yaw rate 
generated by the NASA adaptive algorithm dropped much earlier as compared to the 
aircraft yaw rate. The optimal algorithm generated a yaw rate following the aircraft yaw 
rate for a longer time but also resulted in larger actuator extensions as compared to the 
NASA adaptive algorithm. The washout of the simulator yaw rate generated by the 
optimal algorithm had a larger overshoot than the one generated by the NASA adaptive 
algorithm. The NASA adaptive algoiithm produced a distorted specific force curve with 
an onset greater than the aircraft which decreased too r^idly, while the optimal 
algorithm was closer in shape to the aircraft but with a large amount of overshoot in 
washout. The outputs are shown in Appendix C, Figures C.22, C.23, and C.24. 

8.5. Braking Algorithm 

A y-acceleration pulse with a magnitude of 10 m/sec^ and duration of 10 seconds 
was employed as the input. The optimal algorithm was employed as the washout 
algorithm. The simulator reached its motion limit once and was successfully braked. 
Some large specific force spikes were generated when the simulator was being braked. 
The brake was then released when the washout algorithm ou^ut was driving the 
simulator toward smaller excursions. The outputs are shown in Figure C.25. 
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9. Conclusions 


When the pitch/surge or roll/sway channel was tested, the specific force output of 
the NASA adaptive dgorithm often went in the wrong direction at the beginning of a 
translational acceleration input. The specific force in the wrong direction resulting from 
the cross-coupling from the simulator tilt to the simulator pilot’s head translational 
motion was reduced but not completely eliminated in the first phase of this project. It 
was found that both the magnitude and duration of the specific force in the wrong 
direction wCTe still too large to be ignored when the input contained some frequency 
components near or higher than 0.5 Hz. As discussed in Section 1, negative motion cues 
should always be avoided if possible and when Ac motion cue onset has high importance. 
The specific force in the wrong direction may generate a bad motion sensation which 
happens at the motion onset, therefore it was desired to eliminate it completely. It was 
found that the optimal algorithm successfully handled this problem. The cross-coupling 
was explicitly expressed in the construction of the system equations when the new 
optimal algorithm was developed. Results showed that the specific force in the wrong 
direction was completely eliminated. The optimal algorithm also generated specific force 
outputs with significantly smaller shape distortions than the NASA adaptive algorithm 
did in the pitch/surge and roll/sway test cases. 

In the heave test case, the NASA adaptive algorithm did not generate an obvious 
specific force drop when there was a large drop of the input while the optimal algorithm 
generated one, which was obviously desired. 
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When the pitch channel and the roll channel were tested, the NASA adaptive 
algorithm generated simulator angular velocity outputs with very nice shapes while the 
outputs of the optimal algorithm had some visible distortion. The optimal filters for the 
pitch channel were designed in the pitch/surge channel design. The filters were tuned 
mainly according to test nms with surge inputs. If the filters for the pitch channel are not 
satisfactory, it may be necessary to re-tune those pitch filters without changing the filters 
for the surge channel. This is also true for the filters for the roll chaimel. 

When the yaw channel was tested, the output of the NASA adaptive algorithm 
and the optimal algorithm were differrait If the nonlinear gains for the two algorithms 
are adjusted so diat the two algorithms drive the simulator to the same amount of 
actuator extensions, the simulator angular velocity output of the NASA adaptive 
algorithm will have higher magnitude but be sustained for shorter time than the output of 
the optimal algorithm. It is not clear which output will result in better simulation effects. 

It was found that in most input cases, to generate the specific force at the 
simulator pilot’s head or the simulator angular velocity with the same magnitude, the 
NASA adaptive algorithm and the optimal algorithm would result in about the same 
amount of actuator extensions. The advantages of the optimal algorithm over the NASA 
adaptive algorithm are that the optimal algorithm eliminated the negative motion cues, 
generated outputs with better shapes in many simulation cases, and did not lose some 
desirable motion cues while the NASA adaptive algorithm did in some input cases. 
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A ppendix A. 

Vestibular Rotation Sensation Model Construction and Analysis 

1. Introduction 

The purpose of this study is to find a vestibular semicircular rotation sensation 
model that is most consistent with experimental results and can be justified by theoretical 
analysis, thus providing a reliable approximation to the rotation sensation function. This 
study is based on the reports presented by the many researchers who worked on the 
semicircular canal and rotation sensation function analysis. The model is first 
constructed step by step. The model parameters are then determined and then the model 
is carefully justified. 

2. Model Construction 

The vestibular rotation sensation model has been well studied by several authors. 
Steinhausen [Al] first developed a linear second order model of canal dynamics to 
explain the observed characteristics of vestibular induced eye movements in fish (pike) in 
1931. This model was further refined by the "torsion-pendulum" model of Van Egmond, 
et al. [A2] in 1949, and is later developed from a systems approach by Mayne [A3]. The 
differential equation for this model is 

I9^+c0, +k0, = la (A-1) 

where I = Moment of inertia of the endolymph 

c = Moment of viscous damping of the endol 3 miph 
k = Moment of elastic restoring force of the endolymph 
6e = Angular displacement of the endolymph with respect to the head 
a = Angular acceleration of the head with respect to an inertial axis 
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The following transfer function is obtained from Equation (A-1); 


^e(s) 

a(s) 


1 

2 c k 
s + -s + — 


I 

k 

2 C 

s + — s +1 


(A-2) 


I I k k 
For an overdamped system. Equation (A-2) can be written as 


I 

k 


9 ,( 5 ) ^ 

a(s) (1 +x,s)(l +T25) 


(A-3) 


c I 

wherex, = -’'^2 =~. “^x, > Xj. Schmid, et al. [A4] show that cupula deflection 
iC c 

and the endolymph angular displacement 0e are related by 


, - R A, - 

<l>c =-2 , =-a0. 


h A 


(A-4) 


cp 


where 


R = Central radius of the canal 


Ac = Cross-sectional area of the canal 
Ap = Cross-sectional area of the cupula 
h = Height of the cupula 

and a is a nondimensional number. Substituting (A-4) into (A-3) results in the tr ans fer 
function relating cupula angular deflection to an angular acceleration input: 




a^ 


(A-5) 


a(s) (l+XjSXl+XjS) 

The numerator in Equation (A-5) is equal to the K in the numerator in the transfer 
frmction given by Zacharias [A5], 

Further studies showed that a complete vestibular rotational perception model is 
more complex than the torsion-pendulum model. Young and Oman [A6] formulated an 
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adaptation operator and cascaded it with the torsion-pendulum model to resolve the 
conflicts between the sensed response predicted by the torsion-pendulum model and the 
perceptual response measured in experiments. The addition of adaptation results in the 
following transfer function: 

1 (A-6) 

a(5) l+x^5 (l+t,5Xl+X25) 

Zacharias [A5] noted the following in regard to an additional lead component of 


the form (1 + Xis): 

“In addition to the adaptation just discussed, there appears to be evidence of 
lead sensitivity in vestibular processing of angular velocity irrformation. In studying 
postural reactions to induced body tilt in 1970, Nashner found it necessary to augment 
the torsion pendulum model with a lead term having a 17 msec time constant, in order to 
fit reflex latencies to large amplitude disturbances. As noted by Ormsby, this type of lead 
behavior is not inconsistent with the vestibular nystagmus frequeruy responses reported 
by Benson, in which a high frequency gain rise was noted, consistent with a lead 
(qterator having a 60 msec time constant. Finally, in their investigation of primary 
afferent response of squirrel monkeys to rotational stimuli, Fernandez and Goldberg 
found that the population average frequency response could be best fit with the inclusion 
of a lead term h^ng a 50 msec time constant. ” 

A model representing both the semicircular canal and the peripheral neuron 
transduction dynamics is now established; 




l+x^5 


(A-7) 


a(s) l+x,5 (l+Xj^Xl+XjS) 

Fernandez and Goldberg [A7] determined the vestibular parameters for the squirrel 
monkey by direct measurement of the afferent nerves due to various angular acceleration 
inputs of different am plitudes and frequencies. Their transfer function relates the afferent 
firing rate of the vestibular nerve to the angular acceleration input; 


AFR(5) ^ x,5 l+x^,5 

a(5) l+x„J (1 +x,5X1+''2‘^) 


(A-8) 
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where the sensitivity k is different than the sensitivity K given in Equation (A-7). 

3. Parameter Determination 

Fernandez and Goldberg [A7] detomined the average transfer function 
parameters for the squirrel monkey as 

AFR(5) _ ^ 80s 1+ 0.049s 

a(s) ” 1 + 80s (l + 5.7sXl + 0.003s) 

where the sensitivity k was estimated as 3.44 spikes sec'Vdeg sec'^ for a constant 
acceleration input. Parameters for man are more difficult to measure because direct 
detection of the afferent nerve outputs of the vestibular system cannot be done and 
therefore most experiments were based on subjective responses or nystagmus tests. 

The objectives of several researchers were to determine the values of the 
parameters based on the subjective response of humans. Van Egmond, et al. [A2] 
reported that Xi and t 2 had values of about 10 seconds and 0.1 seconds respectively for 
man in 1949. The values were based on the verbal response of highly trained subjects 
subjected to various motion inputs in both a rotating chair and a torsion swing. Meiry 
(from Zachanas [A5]) measured detection latency as a function of angular acceleration 
step size, and found a long time constant of 7 seconds for roll-axis rotation about the 
earth-vertical axis. Guedry (from Zacharias [A5]) used a short period rotational stimulus 
consisting of an acceleration pulse doublet, and a response measure of apparent 
displacement, and found values of 16 seconds and 7 seconds in yaw and pitch, 
respectively about an earth-vertical axis. Malcolm and Melvill Jones (from Zacharias 
[A5]) investigated the response to earth-vratical rotation about all three body axes by 
using a velocity step as the stimulus, and measured the elapsed time to zero perceptual 
response. They also measured the slow phase velocity (SPV) nystagmus of vestibular 
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induced compensatory eye movements, thus providing two separate measures of canal 
ftmciion. They were able to derive two sets of long time constants for the three axes with 
the torsion pendulum model driving both perceived angular velocity and eye velocity. 
Their results are summarized in Table A. 1 


Table A.I. Torsion Pendulum Long Time Constants (from Zacharias [A5]). 



Yaw 

Pitch 

Roll 

subjective sensation 

10.2 ±1.8 

5.3 ±0.7 

6.1 ±12 

nystagmus SPV 

15.6 + 12 

6,6 ± 0.7 

4.0 ±0.4 


It should be noted that all the above estimations were based on the model expressed as 
the transfer flmction in Equation (A-5). 


Young and Oman [A6] later fit these measured parameters to a more complex 
model. Young and Oman's model also resolved the apparent inconsistency between the 
nystagmus and subjective response measures for yaw rotation. Figure A.1 shows Young 
and Oman's model. The torsion-pendulum model of the canals drives both perceptual 
and nystagmus response channels. Each chaimel has its own adaptation time constant. A 
single long time constant of 16 seconds is used for the torsion-pendulum model. A 
combination of the torsion-pendulum model and the adaptation mechanism generated a 
different response to a step input for each chaimel. When the model is stimulated by step 
inputs, the apparent time constant of the nystagmus decay is 16 seconds, whereas the 
apparent time constant for sensation decay is 10 seconds. This model predicts a result for 
earth-vertical yaw rotation consistent with the data given in Table A.1, but was never 
extended for roll and pitch rotation. 

It can be inferred that the long time constant Xi measured by Van Egmond, et al. 
[A2], Meiry and Malcolm and Melvill-Jones (from Zacharias [AS]) does not actually 
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represent the semicircular canal parameter in the model, but is an overall dynamics 
parameter representing the decay speed of the rotational sensation to a step angular 
velocity input. Zacharias [A5] suggests each axis of rotation has an equivalent “body 
axis” canal pair with a distinct time constant. These “body axis” time constants can then 
be transformed into the three physical canal pairs. The psychophysical results show each 
of the three canal pairs having a distinct value for Xi, especially for yaw. However, 
physiological results based on afferent responses by Fernandez and Goldberg [A7] show 
the same value for ti for all three canal pairs. Zacharias [A5] suggests the differences 
shown in the psychophysical results may occur at a central origin at the perceptual level. 


Nystagmus 120 sec 3.8 deg/sec 

Gain Adaptation Threshold 



0.3 sec Dead Sizbjective 30 sec 1.5 d^sec Scnsaticm of 
Time Delay Gain Adaptation Threshold Vdocity 


Figure A.1, Adaptation Model for Earth-Vertical Rotation (from [A3] ). 
Fernandez and Goldberg [A7] made the following observation in regard to the 
short time constant T 2 : 
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“Our observations, since they extended only to 8.0 Hz. are insufficient to provide 
a direct experimental measure of x 2 . Some estimate of this constant can, however, be 
made from hydrocfynamic considerations. What seems to be required is a solution of the 
Navier-Stokes equation for the complicated geometry represented by the canal, the 
associated ampulla, and the utriculus. No one to our knowledge has accomplished this, 
though Steer has solved the equation for a straight tube in 1967. The cqrproximate value 
of the time constant, so derived, is X 2 = (9 / ‘n where p is the first zero of the 

zeroth-order Bessel function of the first kind; p cmd t) are, respectively, the density and 
viscosity of the endolymph, and r is the internal radius of the tube, in this case the radius 
of the membranous canal. This radius has been measured by Igarashi and leads to a 
value of \2 equal to 0.005 seconds in man and 0.003 seconds in the squirrel monkey . " 

Several authors did experiments to detennine the value of - Their findings are 
summarized in Table A.2. 

Table A.2. Lead Sensitivity to Rotational Stimuli (From Zacharias [A5]). 


(sec) 

Measure 

Source 

0.017 (man) 

posture control 

Nashner 

0.06 (man) 

nystagmus 

Benson and Ormsby 

0.05 (squirrel monkey) 

primary afferent 

Fernandez and Goldberg 


From experiments with human subjects. Van Egmond, et al. [A2] showed that the 
perceived angular velocity © is proportional to the cupula deflection <j>c by the long time 
constant xi, resulting in the transfer function 


^ = (A-IO) 

©( 5 ) (l+X,j)(l+t25) 

As shown in Equation (A-9), Fernandez and Goldberg [A7] show a gain sensitivity k 
between the input stimulus and the afferent firing rate that was estimated at 3.44 spikes 
sec'* / deg sec'^. Ormsby (from Zacharias [A5]) proposed that the perceived angular 
velocity ® is proportional to the afferent firing rate. While no one to date has 
experimentally obtained this parameter, Zacharias [A5] noted that Curry, et al. provided 
an estimate of the overall gain between perceived and input angular velocity based on 
angular acceleration thresholds. 
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From subjective pilot measurements of angular acceleration thresholds on a 
moving base platform, Hosman and van der Vaart [A8] show the perceived input 
threshold amin as 


amin = 


^min 

|h(<d)| 


(A-ll) 


where bniin = Minimum sensed amplitude of the afferent firing rate 

1h(©)| = Modulus of the vestibular transfer function 
From the measured frequency response of a^n the following transfer function is obtained, 
neglecting gain sensitivity and adaptation: 


H(s) = 


1 + 0.1097 5 


(1 + 5.924 5)0+0.0055) 


(A-12) 


These results are based upon roll and pitch acceleration thresholds; yaw thresholds were 
not measured. The value for ti agrees well with the value obtained by Fernandez and 
Goldberg. The value obtained for tl is nearly twice the nystagmus value obtained by 
Benson and Ormsby as given in Table A.2. 

From the results reported by the authors mentioned in this section, a transfer 
function that can best describe the vestibular rotational sensation system is proposed; 

AFR(5) 805 1 + 0.065 

— = 3.44 (A-13) 

a(5) 1 + 805 (1 + 5.735)(1 + 0.0055) ^ ^ 

It should be noted that the parameters in the above transfer function are by no means 

exact, but are of the correct order of magnitude, thus making the transfer function a 

meaningful approximation to the real d 3 mamics of the vestibular rotation sensation 

system. The frequency response of the transfer function given in Equation (A-13) with 
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gain k = 1 is shown in Figure A.2. Both the torsion-pendulum model and the complete 
sensation model with lead and adaptation mechanisms included are both shown. 


Bode Plot of Rotational Sensation TF w/Ang Accel Input 



Frequency (Hz) 


Figure A^. Frequency Response of Vestibular Rotation Sensation System. 

The sensory function of the semicircular canal can be described by observing the 
frequency response of the torsion-pendulum model. In the range of normal head 
movement from 0.05 to 5.0 Hz (Mayne, [A3]), the gain response decreases by 20 
dB/decade with the phase close to minus 90 degrees. In this frequency range the canal 
functions as an “integrating accelerometer” or an angular velocity transducer. At very 
low frequencies less than 0.01 Hz, the phase approaches zero degrees, thus functioning as 
an accelerometer. At very high frequencies greater than 100 Hz, the phase approaches 
minus 180 degrees, thus functioning as an angular displacement transducer. The effects 
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of adaptation and lead on rotational sensation are apparent; adaptation influences low 
frequencies below 0.01 Hz while lead influences high frequencies greater than 10 Hz. 

4. Physiological Interpretation 

Modem theories of the operation of the semicircular canal receptors are based on 
the assumption that the nervous impulses are generated by deflection of the hairs in the 
sensory cells as a result of the cupula displacement. Input acceleration is first 
transformed into cupula deflection by the cupula-endolymph system. Then the cupula 
deflection is further transformed into electrical impulses by the mechano-neural 
transduction system consisting of sensory hair cells, afferent nerves and efferent nerves. 

The cupula-endolymph system was likened to an overdamped linear torsion- 
pendulum by Steinhausen [Al] in 1931. A rigorous analytical evaluation of the dynamics 
of the endolymph motion in the semicircular canals was made by Steer [A9] in 1967. 
The torsion-pendulum model is well accepted by different authors without controversy 
upon the form of the model. 

Besides the torsion-pendulum dynamics there are some additional terms in the 
complete transfer flmction in Equation (A-8). These terms could be grouped as 


kT,5(l+T^5) 

O+'c.^) 


(A-14) 


representing an adaptation-lead mechanism. Controversy occurred when different 
authors tried to interpret this adaptation-lead mechanism primarily because the precise 
mechanism of hair cell stimulation is not yet understood. 

The first controversy is whether the adaptation arises in the mechanics of the 
cupula-endolymph system or in the mechano-nmral system. Goldbra^g and Fernandez 
[A8] presented a good discussion of the origin of the adaptation: 
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“Two observations tend to suggest that the adcptation arises in the mechano- 
neural system. The first is the clear adaptation seen by Lowenstein when polarizing 
currents were applied to the vestibular nerve of the thomback ray. Presumably, the 
currents acted Erectly on the nerve terminals. The second is the fact that units differ 
greatly in their adaptive properties. Were adaptation of mechanical origin, one would 
have to assume that hair cells differ in the way they are mechanically coupled to the 
motion of the cupula An assumption which appears to us more reasonable is that the 
adcptation reflects the physiology of the hair cells catd/or of the efferent nerve terminals. 
Another, perheps unlikely, possibility is that the adeptation results from the activation of 
the efferents innervating the sensory epithelium. " 

Controversy arose again about the origin of the lead operator. Since no 
mechanism was found in the cupula-endolymph system that could provide a reasonable 
interpretation for the lead operator, most authors suggested that the lead operator arose 
from the mechano-neural system. Fernandez and Goldberg [AT] su gg ested that the lead 
operator implies the sensory hair cells are sensitive to both the displacement and the 
velocity of the cupula. The time constant tl reflects the relative sensitivities to these two 
aspects of the cupular motion, 
mechanism represents the transfer 
cupula deflection ()>; 

F(^) 

(j)(5) 

This expression implies that both 
system. 

Schmid, Buizza, et al. [A4] gave two other plausible interpretations for the 
adaptation-lead mechanism. One is based on the observation that there are two types of 
sensory cells in the vestibular sensory epithelia of mammals. All the sensory cells are 
assumed to be characterized by a transfer flmetion of the type 


They flirtho^ suggested that the adaptation-lead 
function between the discharge frequency F and the 


1cT,jO+^l) 

1 + 

adaptation and lead arise from the mechano-neural 
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(A-16) 


F,(s) K*T*s 
(^(s) 1 + T*s 


where F, is the afferent nerve output, (|) is the cupula deflection, Ki is the sensitivity 
factor and T* is the time constant of the transduction process. It is further assumed that 
the two different types of sensory cells work in parallel with a different T*. Then the 
overall transfer function of the mechano-neural transduction system can be expressed as 


F(s)_ K„’i;s ^ K^T,s 
(|)(s) 1 + T,s 1 + TjS 

By simple manipulation the following transfer function can be obtained; 

F(s)_ K^JisQ + T^s) 
m (1 + T,s)(l + T,s) 


(A-17) 


(A-18) 


where T, » and Tj » 2 X 3 . If T 2 is assumed negligible, then the transfer function can 
be reduced to an adaptation-lead operator and give a good interpretation to the origin of 
the adaptation-lead mechanism. A very small Tj implies that a group of sensory cells 
have a very small adaptation time constant. Experiments have shown that the adaptation 
time constants of the vestibular mechano-neural system of squirrel monkeys range 
approximately from 30 seconds to infinitely long. The fact that no group of sensory cells 
with a small adaptation time constant were found exp^imentally conflicts with the first 
interpretation presented by Schmid, et al. [A4]. 

The second interpretation presented by Schmid, et al. [A4] is based on the 
existence of efferent pathways descending from the vestibular nuclei to the sensory 
epithelium, possibly with an inhibitory function. These efferents may represent a 
negative feedback that modifies the characteristics of the mechano-neural transduction. 
The average dynamics of all sensory cells is assumed to be 
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(A-19) 


F,(s) K*T*s 
(|)(s) ~ 1 + T*s 


The dynamics of the inhibitory process is described by means of a first order system with 
a time constant T. much smaller than T * . The complete mechano-neural transduction 
system can be represented by the block diagram in Figure A.2. 


m 




^ K*T*s 
\ 1 + T*s 

Dynamics of the 
sensory cells 


U(s) 


K. 


1 + Ts 


Dynamics of the 
efferent pathway 

Figure A.2. A Closed-loop Interpretation of the Mechano-neural System. 
The closed-loop transfer function can be obtained: 


K*T*s 
1 + T*s 


F(s) 

<Ks)-^, K*T*s K. 


K*T*s(l + TjS) 


1 + T*s 1 + Ts 


1 + (K * K; + 1)T * s + TjS + T * TjS" 


(A-20) 


By introducing the simplifications suggested by the condition T^ « T*, it follows that 


F(s) 


K*T*s (1+T^s) 


[l+(K*Ki+l)T*s](l+j^;^s) 


(A-2I) 


If 


is small enough to be considered negligible, then the above transfer 


K*K.+1 

function can be further reduced to the form of Equation (A- 15). 
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4. Conclusion 


The vestibular rotational sensation system has been well studied. Its model 
consists of a second order torsion-pendulum operator and an adaptation-lead operator. 
The torsion-pendulum operator arises from the cupula-endolymph system and functions 
as an angular velocity transducer in the range of normal head movement. The model is 
then augmented by the adaptation-lead operator that arises from the mechano-neural 
transduction system. The adaptation-lead operator most likely represents a closed-loop 
dynamic system consisting of sensory cells, afferent nerves as feed-forward path and 
efferent nerves as negative feedback path. Based upon the model parameters reported in 
the literature a transfer function best approximating the rotational sensation system is 
proposed. 
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Appendix B. 


Phase 1 Rotational 

Output Figures 

Figure B.l. 

Classical Algorithm Roll Doublet Pulse Input. 
1 m/s^ peak, 10 second duration input. 

Figure B.2. 

Optimal Algorithm RoU Doublet Pulse Input. 
1 m/s^ peak, 10 second duration input. 

Figure B.3. 

NASA Adaptive Algorithm Roll Doublet Pulse Input. 
1 m/s^ peak, 10 second duration input. 

Figure B.4. 

UTIAS Adaptive Algorithm RoU Doublet Pulse Input. 
1 m/s^ peak, 10 second duration input. 
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Figure B.l. Classical Algorithm Roll Doublet Pulse Input. 
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ASF-Z & SSF-Z 



Figure B.I. Classical Algorithm Roll Doublet Pulse Input. 
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Figure B.l. Classical Algorithm Roll Doublet Pulse Input. 
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Figure B.2. Optimal Algorithm Roll Doublet Pulse Input. 
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Figure B.2. Optimal Algorithm Roll Doublet Pulse Input. 
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Figure B.2. Optimal Algorithm Roll Doublet Pulse Input. 
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Figure B.2. Optimal Algorithm Roll Doublet Pulse Input. 
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Figure B.3. NASA Adaptive Algorithm Roll Doublet Pulse Input. 
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Figure B3. NASA Adaptive Algorithm Roll Doublet Pulse Input. 
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Figure B.3. NASA Adaptive Algorithm Roll Doublet Pulse Input. 
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Figure B^. NASA Adaptive Algorithm Roll Doublet Pulse Input. 
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Figure B3. NASA Adaptive Algorithm Roll Doublet Pulse Input. 
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Figure B.4. UTIAS Adaptive Algorithm Roll Doublet Pulse Input. 
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Figure B.4. UTIAS Adaptive Algorithm Roll Doublet Pulse Input. 
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Figure B.4. UTIAS Adaptive Algorithm Roll Doublet Pulse Input. 
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Figure B.4. UTIAS Adaptive Algorithm Roll Doublet Pulse Input. 
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Figure B.4. UTIAS Adaptive Algorithm Roll Doublet Pulse Input. 
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Figure B.4. UTIAS Adaptive Algorithm Roll Doublet Pulse Input. 
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Figure B.4. UTIAS Adaptive Algorithm Roll Doublet Pulse Input. 
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A ppendix C. 

Phase 2 Test Run Output Figures 


Figure C.l. 

Optimal Algorithm Surge Half Sinusoidal Pulse Input. 
0. 1 Hz half cycle, 1 m/s^ peak x-accderation input. 

Figure C.2. 

NASA Adaptive Algorithm Surge Half Sinusoidal Pulse Input. 
0. 1 Hz half cycle, 1 m/s^ peak x-acceleration iiq)ut. 

Figure C.3. 

Surge Pulse Input Specific Force Comparison. 
Comparison of the specific forces for both algorithms. 

Figure C.4. 

Optimal Algorithm Surge 3 m/sVs Ramp-Step Input. 
3 m/s^/s slope, 5 m/s^ peak x-acceleration input. 

Figure C.5. 

NASA Adaptive Algorithm Surge 3 m/s^/s Ran^Step Ii^ut. 
3 m/s^/s slope, 5 m/s^ peak x-accderation input. 

lagure C.6. 

Surge 3 m/s^/s Ramp-Step Input Specific Force Comparison. 
Conq)arison of the specific forces for both algorithms. 

Figure C.l. 

Optimal Algorithm Surge 5 m/sVs Ramp-Stq> Input. 
S m/s^/s slope, S m/s^ peak x-acceleration input. 

Figure C.8. 

NASA Adaptive Algorithm Surge 5 m/s^/s Ramp-Step Irq)ut. 
5 m/s^/s slope, 5 m/s^ peak x-accderation irq>ut. 

Figure C.9. 

Surge 5 m/sVs Ramp-Step Irqnit Specific Force Comparison. 
Comparison of the specific forces for both algorhhms. 

Figure CIO. 

Optimal Algorithm Sway Half Sinusoidal Pulse Ir^ut. 
0.5 Hz half cycle, 1 m/s^ peak y-acceleration iiq)ut. 

FigureC.il. 

NASA Adaptive Algorithm Sway Half Sinusoidal Pulse Input. 
0.5 Hz half cycle, 1 m/s^ peak y-acceleration input. 

Figure C. 12. 

Sway Pulse Input Specific Force Comparison. 
Comparison of the spedfic forces for both algorithms. 
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Figure C. 13. 

Optimal Algorithm Heave Square Pulse bq)ut. 

3 m/s^ peak, 5 second duration z-acceleration iiqiut. 

Figure C. 14. 

NASA Adaptive Algorithm Heave Square Pulse Input. 
3 m/s^ peak, 5 second duration z-acceleration input. 

Figure C.15. 

Heave Square Pulse Input Spedfic Force Conq)aiison. 
Comparison of the specific forces for both algorithms. 

HgureC.16. 

Optimal Algorithm Roll Doublet Pulse Iiqnit. 
0. 1 rad/sec^peak, 5 second duration input. 

Figure C. 17. 

NASA Ad^tive Algorithm Roll Doublet Pulse Input. 
0. 1 rad/sec^ peak, 5 second duration iiqiut. 

HgureC.18. 

Roll Doublet Spedfic Force Comparison. 

Comparison of the spedfic forces for both algorithms. 

RgureC.19. 

Optimal Algorithm Pitch Doublet Pulse Iiqmt 
0. 1 rad/sec^ peak, 5 second duration input. 

Figure C.20. 

NASA Adaptive Algorithm Pitch Doublet Pulse Input. 
0.1 rad/sec^ peak, 5 second duration iiq}ut 

Rgure C.21. 

Pitch Doublet Spedfic Force Cooqiarison. 
Comparison of the specific forces for both algorithms. 

Figure C.22. 

Optimal Algorithm Yaw Doublet Pulse Input. 
0. 1 rad/sec^ peak, 5 second duration input. 

I^gure C.23. 

NASA Adaptive Algorithm Yaw Doublet Pulse Input 
0. 1 rad/sec^ peak, 5 second duration iiqiut. 

Hgure C.24. 

Yaw Doublet Spedfic Force Comparison. 

Comparison of the spedfic forces for both algorithms. 

Hgure C.25. 

Sway Pulse Input Braking Algorithm Test. 

10 m/s^ peak, 10 second duration y-accderation input. 


Spedfic force comparison figures are labelled as follows: 

Line with no marks — aircraft response. 

Line with marics — simulator response: o — NASA Adaptive; * — Optimal 
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Figure C.l. Optimal Algorithm Surge Half Sinusoidal Pulse Input. 


153 








PsI (deg) Theta (deg) Phi (deg) 



t(sec) 


t(sec) 


Figure C.l. Optimal Algorithm Surge Half Sinusoidal Pulse Input. 
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Figure C.l. Optimal Algorithm Surge Half Sinusoidal Pulse Input. 
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Figure C.l. Optimal Algorithm Surge Half Sinusoidal Pulse Input. 
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Figure C.1. Optimal Algoiithm Surge Half Sinusoidal Pulse Input. 
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Figure C.l. Optimal Algorithm Surge Half Sinusoidal Pulse Input. 
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Figure C.2. NASA Adaptive Algorithm Surge Half Sinusoidal Pulse Input. 
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Figure C^. NASA Adaptive Algorithm Surge Half Sinusoidal Pulse Input. 
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Figure C^. NASA Adaptive Algorithm Surge Half Sinusoidal Pulse Input. 
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Figure C.2. NASA Adaptive Algorithm Surge Half Sinusoidal Pulse Input. 
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Figure NASA Adaptive Algorithm Surge Half Sinusoidal Pulse Input. 
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Fi^re CJ2. NASA Adaptive Algorithm Surge Half Sinusoidal Pulse Input. 
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Figure C.4. Optimal Algoritiim Surge 3 m/s^/s Ramp-Step Input. 










AW-r (deg/s) AW-q (deg/s) AW-p (deg/s) 


Aircraft Ar^ular Rate Simulator Angular Rate 




Figure C.4. Optinial Algorhfam Surge 3 m/sVs Ramp-Step Input. 
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Figure C.4. Optimal Algorithm Surge 3 m/sVs Ramp-Step Input. 


171 










t (sec) t (sec) 

Figure C.4. Optimal Algorithm Surge 3 m/s^/s Ramp-Step Input. 
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Figure C.4, Optimal Algorithm Surge 3 m/sVs Ramp-Step Input. 
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Figure C^. NASA Adaptive Algorithm Suige 3 m/sVs Ramp-Step Input. 
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Figure C^. NASA Ad^rtive Algorithm 1 
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Figure C^. NASA Ad^)tive Algoiithm Surge 3 m/s^/s Ramp-Step Input. 
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Figure C5. NASA Adaptive Algorithm Surge 3 m/s^/s Ramp-Step Liput. 
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Figure C5. NASA Adqrtive Algorithm Surge 3 m/s^/s Ramp-Step Input. 
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Figure C.6. Surge 3 m/sVs Ramp-Step Input Specific Force Comparison. 
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Figure C.7. Optimal Algorithm Surge 5 m/s^/s Ramp-Step Input. 
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Figure C7. Optimal Algorithm Surge 5 m/^/s Ramp-Stq? Input. 
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Figure C7. Optimal Algorithm Surge 5 m/s^/s Ramp>Step Input. 
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Figure C.7. Optimal Algorithm Surge 5 m/s^/s Ramp-Step Input. 
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figure C.7. Optimal Algorithm Surge 5 m/sVs Ramp-Step Input. 
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Figure C8. NASA Ad^rtive Algorithm Surge 5 m/s^/s Ramp-Step Input. 
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Figure C8. NASA Adaptive Algorithm Surge 5 m/sVs Ranq)-Step Ii^ut 
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Figure C^. NASA Adaptive Algorithm Surge S m/sVs Ranq)-Step Input. 


194 










Simu. Translertional Displacement 









Jack 5 (m) Jack 3 (m) Jack1(m) 


Jack Extension 


Jack Extension 





t(sec) 





Figure C.8, NASA Adaptive Algorithm Suige 5 m/s^/s Ramp-Step Ir^ut. 
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Figure C.10. Optimal Algorhbm Sway Half Sinusoidal Pulse Input. 
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Figure CIO. Optimal Algorithm Sway Half Sinusoidal Pulse Input. 
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Figure C.10. Optimal Algorithm Sway Half Sinusoidal Pulse Input. 
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Figure C.10. Optimal Algorithm Sway Half Sinusoidal Pulse Input. 
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Figure C.11. NASA Adaptive Algorithm Sway Half Sinusoidal Pulse Input. 
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Figure C.11. NASA Adaptive Algorithm Sway Half Sinusoidal Pulse Input. 
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Figure C.11. NASA Additive Algorithm Sway Half Sinusoidal Pulse Input. 
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Figure C.11. NASA Adaptive Algorithm Sway Half Sinusoidal Pulse Iiq)ut. 
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Figure C.11. NASA Ad^tive Algorithni Sway Half Sinusoidal Pulse Input. 
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Figure C.12. Sway Half Sinusoidal Pulse Input Specific Force Comparison. 
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Figure C.13. Optimal Algorithm Heave Square Pulse Iiq)ut. 
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Figure C.13. Optimal Algonthm Heave Square Pulse Input. 
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Figure C.13. Optimal Algorithm Heave Square Pulse Input 
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Figure C.13. Optimal Algorithm Heave Square Pulse Input. 
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Figure C.13. Optimal Algorithm Heave Square Pulse Liput. 
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Figure C.14. NASA Adaptive Algoiithin Heave Square Pulse Input. 
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Figure C.14. NASA Adaptive Algorithm Heave Square Pulse Input. 
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Figure C.14. NASA Adaptive Algorithin Heave Square Pulse Ii]q)ut. 
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Figure C. 14 . NASA Adaptive Algorithm Heave Square Pulse Iiq)ut. 
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Figure C.14. NASA Adaptive Algoiithm Heave Square Pulse Input. 
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Figure C.15. Heave Square Pulse Input Specific Force Comparison. 
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Figure C.15. Heave Square Pulse Input Specific Force Comparison. 
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Figure C.16. Optimal Algorithm Roll Doublet Pulse Input. 
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Figure C.16. Optimal Algorithm Roll Doublet Pulse Input. 
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Figure C.16. Optimal Algoritfam Roll Doublet Pulse Input. 
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Figure C.16. Optimal Algorithm Roll Doublet Pulse Input. 
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Figure C.16. Optimal Algorithm Roll Doublet Pulse Input. 


233 















Jack S (m) Jack 3 (m) Jack 1 (m) 



t(sec) 


t(sec) 


Figure C.16. Op timal Algorithm Roll Doublet Pulse Input. 
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Figure C17. NASA Adaptive Algorithm Roll Doublet Pulse Input. 
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Figure C.17. NASA Adaptive Algorithm Roll Doublet Pulse Input. 
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Figure C.17. NASA Adaptive Algorithm Roll Doublet Pulse Input 
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Figure C.18. Roll Doublet Pulse Input Spedfic Force Comparison. 
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Figure C.19. Optimal Algorithm Pitch Doublet Pulse Input. 
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F^ure C.19. Optimal Algorithm Pitch Doublet Pulse Input. 
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F^re C.19. Optimal Algoritim Pitch Doublet Pulse Input. 
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Figure C20. NASA Adaptive Algorithm Pitch Doublet Pulse Input. 
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Figure C20. NASA Adaptive Algorithm Pitch Doublet Pulse Input. 
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Figure C20. NASA Adaptive Algorithm Pitch Doublet Pulse Input. 
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Figure C^O. NASA Additive Algorithm Pitch Doublet Pulse L]^ut. 
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Figure C^O. NASA Ad^)tive Algoiithni Pitch Doublet Pulse Input. 
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Figure C21. Pitch Doublet Pulse Input Spedfic Force Comparison. 
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Figure C22. Optimal Algorithm Yaw Doublet Pulse Input. 
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Figure C22. Optinial Algorithm Yaw Doublet Pulse Input. 
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Figure C22. Optimal Algorithm Yaw Doublet Pulse Input. 
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Figure CJ3. NASA Adaptive Algoiitlim Yaw Doublet Pulse Iiq)ut. 
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F^re CJ23, NASA Adaptive Algorithm Yaw Doublet Pulse Input. 
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Figure C^. NASA Adaptive Algorithm Yaw Doublet Pulse Input. 
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Figure C25. Sway Pulse Input Braking Algorithm Test. 
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F^re C25. Sway Pulse Input Braking Algorithm Test. 
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Figure C25. Sway Pulse Input Braking Algorithm Test. 
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Figure C25. Sway Pulse Input Braking Algorithm Test. 
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Figure C25. Sway Pulse Input Braking Algorithm Test. 
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